République Algérienne Démocratique et Populaire 
J <>dl fokUl SjIJj 

Ministère de l’enseignement Supérieur et de la Recherche scientifique 



Université Mohamed Khider Biskra 
Faculté des Sciences et de la Technologie 
Département de Génie Electrique 
Filière : Electronique 

Option : Microélectronique et hyperfréquences 


Réf : 


Mémoire de Fin d' Etudes 
En vue de l’obtention du diplôme : 

MASTER 

Thème 


J 


Etude par simulation numérique 
d'une cellule solaire en CIGS 




Présenté par : 

Kalache Salim 

Soutenu le : 01 Juin 2017 

Devant le jury composé de : 

M me Terghini Warda M.C.B 

M me Tobbeche Souad Pr. 

M me Ouarhlent Saloua M.A.A 


Président 

Encadreur 

Examinateur 




Année universitaire : 2016 / 2017 




République Algérienne Démocratique et Populaire 
^bül Sjljj 

Ministère de l’enseignement Supérieur et de la recherche scientifique 



Université Mohamed Khider Biskra 
Faculté des Sciences et de la Technologie 
Département de Génie Electrique 
Filière : Electronique 

Option : Microélectronique et hyperfréquences 


Mémoire de Fin d'Etudes 
En vue de l’obtention du diplôme : 

MASTER 

Thème 

Etude par simulation numérique d'une cellule 

solaire en CIGS 


Présenté par : Avis favorable de l’encadreur : 

Kalache Salim Tobbeche Souad Signature 


Avis favorable du Président du Jury 

Terghini Warda 


Cachet et signature 


République Algérienne Démocratique et Populaire 
diaJlj <^bdl jaakjil Sjljj 

Ministère de l’enseignement Supérieur et de la Recherche scientifique 



Université Mohamed Khider Biskra 
Faculté des Sciences et de la Technologie 
Département de Génie Electrique 
Filière : Electronique 

Option : Microélectronique et hyperfréquences 

Thème : 


Etude par simulation numérique d'une cellule 

solaire en CIGS 


Proposé par : Pr. Tobbeche Souad 
Dirigé par : Pr. Tobbeche Souad 


Résumé 


Le travail de ce mémoire de master est une étude par simulation numérique des caractéristiques électriques 
d’une cellule solaire en couches minces et à hétéroj onction à base de CIGS. En premier lieu, nous avons 
simulé deux cellules solaires de structures : ZnO/CdS(n)/CIGS(p)/Mo et ZnO/ZnS (n)/CIGS(p)/Mo en 
utilisant le logiciel Silvaco Atlas-2D. La cellule solaire avec la couche tampon ZnS donne le meilleur 
rendement i\=23.54% par rapport à la cellule solaire avec la couche tampon en CdS, i\=22.60%. En 
deuxième lieu, nous avons optimisé la cellule solaire en ZnS/CIGS en étudiant l’effet du dopage et de 
l’épaisseur des couches en ZnS et CIGS. Le rendement optimal obtenu est de 26.82% pour des épaisseurs 
des couches en ZnS et CIGS de l’ordre de 25 nm et 4 pm respectivement et dopés par des concentrations de 
l’ordre 6 x 10 17 et 10 18 cm" 3 respectivement. 
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Introduction générale 


L’énergie est l’un des moteurs de développement des sociétés, elle est aussi, le pilier de 
l’économie moderne. Elle est l’étincelle qui permet l’accomplissement de toute activité 
humaine. Ses sources se sont diversifiées au cours du temps afin de satisfaire les besoins 
toujours accrus de l’industrie et des consommateurs, une des solutions les plus prometteuses 
pour la future énergétique de l’humanité, est le photovoltaïque (cellule solaire). 

La cellule solaire est un composant électronique qui peut convertir la lumière du soleil en 
électricité, cette conversion se produit au sein de matériaux semi-conducteurs, qui ont comme 
propriété de libérer leurs porteurs de charge (électrons et trous) sous l’influence d’une 
excitation extérieure. Actuellement l’objectif de multiples recherches est de réaliser la 
meilleure cellule avec un haut rendement énergétique et un prix de revient bas. Pour diminuer 
à la fois le poids et le prix de ces dispositifs tout en augmentant leurs surfaces et leurs 
souplesses, des cellules en couches minces ont été réalisées. 

Les cellules solaires en couches minces qui utilisent le séléniure de Cuivre, d’indium et de 
Gallium (CIGS) représentent l’approche la plus prometteuse au plan de réduction des coûts de 
production. L’avantage de ce matériau est son élaboration facile d’une part et excellente 
stabilité et une grande résistance aux rayonnements d’autre part [1, 2], cette cellule est 
caractérisée actuellement par un rendement très élevé de 22.6% et sa durée de vie est très élevée 
[3-5], 

Le sulfure de cadmium (CdS) est largement utilisé comme couche tampon pour les cellules 
solaires en couche mince à base du CIGS. Le gap optique de CdS est de 2.42 eV, ce qui offre 
une large fenêtre d'absorption du rayonnement solaire sur la couche d'absorption en CIGS. Plus 
récemment, la recherche sur les cellules solaires en CIGS sans utiliser le matériau CdS a été 
intensifiée pour réduire les coûts de production et de recyclage et pour réduire la toxicité du 
cadium (Cd) [6, 7], Compte tenu de ces facteurs, il y’ a une compatibilité de la couche 
d'absorption en CIGS avec d'autres couches tampons à large bande interdite. Le sulfure de zinc 
(ZnS) préparé par dépôt de bain chimique (CBD) s'avère être une alternative attrayante au CdS 
en combinaison avec des absorbeurs de CIGS [8, 9] pour leur bande interdite large de l’ordre 
de 3.68 eV et non toxique pour l’environnement. 
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Introduction générale 


A cause des coûts élevés de l’expérimentation, les chercheurs se sont dirigés à l’utilisation 
de logiciels de simulation numérique telle que SILVACO pour prédire la réponse des structures 
à base de matériaux préparés pour la conversion photovoltaïque. Ces logiciels permettent des 
centaines de combinaisons de paramètres avant la fabrication réelle du dispositif. C’est un 
moyen rapide, consistant et pas cher pour le design de cellules solaires. Il permet, entre autre, 
la visualisation des paramètres qui influencent le comportement des cellules solaires. 

L’objectif de ce présent mémoire est l’étude par simulation numérique de la cellule solaire 
à base de Cu(In,Ga)Se2 pour deux couches tampon en CdS et en ZnS. Nous allons aussi simuler 
l'effet de l'épaisseur et le dopage des couches en ZnS et CIGS en vue d’optimiser la cellule 
solaire en CIGS. 

Pour ce faire, notre travail est divisé en quatre chapitres comme suit : 

Z Dans le premier chapitre : nous allons présenter les propriétés physiques des 
semiconducteurs en générales et le semi-conducteur CIGS en particulier constituant 
la cellule solaire, et procédés d’élaboration de ses cellules. 

Z Le deuxième chapitre est une étude sur les cellules solaires en général et les cellules 
solaires en CIGS en particulier. Nous décrirons le principe de fonctionnement et les 
caractéristiques de ces cellules solaires. 

Z Le troisième chapitre présente le logiciel Silvaco utilisé pour la simulation de la cellule 
solaire. 

Z Et enfin, Dans le dernier chapitre nous présenterons et discuterons les résultats 
obtenus. 
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Chapitre I 


Propriétés Du Matériau CIGS 


1.1 Introduction 

Beaucoup de laboratoires de recherche à travers le monde travaillent sur les semi- 
conducteurs chalcopyrites pour comprendre les propriétés électroniques, structurales, 
élastiques, mécaniques et optiques de ces matériaux. Cet engouement à ces matériaux est dû à 
leurs propriétés optiques et électroniques intéressantes et à leurs diverses applications pratiques 
dans les cellules solaires, les détecteurs photovoltaïques et l’électronique. Ils servent de base 
surtout pour la cellule solaire ainsi que les modules et capteurs photovoltaïques. La recherche 
dans ce domaine a toujours comme but majeur d’améliorer le rendement électrique et quantique 
de ces cellules photo voltaïques à base de couches minces. 

Le CIGS dans sa structure chalcopyrite possède des propriétés électroniques, électriques et 
optiques intéressantes. La bande interdite est de nature directe et son coefficient d’absorption 
qui est très élevé et de l’ordre de 10 5 cm' 1 dans la gamme du spectre solaire. Le composé CIGS 
est considéré comme matériau semiconducteur important pour les applications 
optoélectroniques et photovoltaïques. 

Dans ce chapitre, nous allons présenter la différence entre conducteurs, isolants et semi- 
conducteurs, les propriétés physiques des semiconducteurs en CIGS constituant la cellule 
solaire et les procédés d’élaboration. 

1.2 Conducteurs - Isolants - Semi-conducteurs 

Les matériaux ayant la plus faible résistivité à température ambiante, typiquement inférieure 
à 10" 5 Qcm, sont les métaux (cuivre, or, argent, aluminium...). La conduction électrique 
s’effectue essentiellement par les électrons libres dont la concentration diffère peu d’un métal 
à l’autre (de 10 22 à 10 23 cm" 3 ) quelle que soit sa pureté. Une augmentation de la température 
provoque une légère augmentation de la résistivité, pouvant s’expliquer par le fait que les 
électrons libres sont gênés dans leur déplacement par les vibrations (croissantes avec la 
température) des atomes du métal. 

Les matériaux dont la résistivité est typiquement supérieure à 10 8 Qcm sont considérés 
comme isolants ; c’est le cas pour le verre, le mica, la silice ( Si0 2 ), le carbone (diamant). Cette 
fois l’augmentation de la température peut provoquer la libération d’électrons (ainsi que de 
"trous") qui peuvent participer à la conduction électrique, ce qui provoque une baisse de la 
résistivité avec la température. 
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Entre les métaux et les isolants se trouvent les semi-conducteurs (SC) dont la résistivité varie 
de 10~ 3 à 10 4 Qcm (ou plus). La conduction électrique se fait par les électrons et les trous, ou 
de façon préférentielle par l’un ou l’autre type de porteurs. Un semi-conducteur peut être soit 
pur auquel cas il est dit "intrinsèque", soit dopé par des impuretés (qui permettent de contrôler 
sa résistivité) auquel cas il est dit "extrinsèque". Si on prend, par exemple, du Silicium assez 
pur et qu’on lui ajoute un atome de Bore ou de Phosphore pour 10 5 atomes de Silicium, sa 
résistivité passe de 10 3 à environ 10 2 Qcm. 

Le tableau 1.1 donne des exemples de matériaux ou de composés semi-conducteurs en 
fonction des éléments qui les constituent et de la position de ces éléments dans le tableau de 
Mendeleïev. 



Tableau 1.1 : Exemples de semi-conducteurs. 
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1.3 Semi-conducteurs intrinsèques 

Un semi-conducteur est dit intrinsèque lorsque le cristal n’est pas pollué (volontairement ou 
non) par des impuretés pouvant changer la concentration en porteurs libres. A température 
constante, un équilibre s’établit (Fig.I.l) entre les phénomènes d’ionisation thermique et de 
recombinaison ; les électrons libres et les ions de silicium apparaissant en quantités égales. La 
concentration en électrons libres n et en trous libres p sont égales à ni la concentration 
intrinsèque. 

La mécanique statistique montre que la population des porteurs libres (n la concentration 
des électrons en cm' 3 dans la bande de conduction et p la concentration des trous en cm' 3 dans 
la bande de valence) s’exprime selon les lois : 

n = N c exp (- ^r) p = N v exp (- 1. 1 

Où N c et N v sont respectivement la densité effective d’états des électrons dans la bande de 
conduction et la densité effective d’états des trous dans la bande de valence. Ces deux 
coefficients évoluent avec la température selon une loi en T 3/2 . 

• A E n et A E p représentent deux différences d’énergies liées à un niveau dit de Fermi qui 

indique les écarts de population entre les électrons et les trous. 


Oj 

'«b 

ï— 
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Bande de conduction 


AE 

n = N r exp(— - ) 

c t' KT 


E C 


E V 


bande interdite 
1.12 eV pour Si 


Bande de valence 



Efî : niveau de Fermi 


\ p = expf — ) 

v KT 


Fig.I.l : populations des porteurs du S.C. intrinsèque et niveau de Fermi. 
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La concentration intrinsèque n* en électrons libres et en trous libres par cm' 3 dépend de la 
hauteur de bande interdite E g et de la température T selon la loi : 


3 

n = p = m = ATz. exp 



1.2 


A : constante du matériau 

Eg : hauteur de bande interdite (eV) 

K : constante de Boltzmann = 8, 6 .10" 5 eV K" 1 
T : température absolue en K 

1.4 Semi-conducteurs extrinsèques 

1.4.1 Semi-conducteurs du type n 

On obtient un S.C. de type n en injectant dans le cristal de silicium des atomes qui possèdent 
5 électrons sur leur couche périphérique (phosphore ou arsenic de la 5° colonne de la 
classification) (figure 1.2). 



Bande de conduction 


r = N d 


T """ T AE n 

Eg — Ejj — L-- 


Ep n : niveau de Fermi 

A E n = ÂYln(— ) 
n- 


Ey! 

Bande de valence 


v_»? 


p= T d 


Fig.I.2 : libération d’un électron par l’atome de phosphore et schéma des bandes. 

A la température ordinaire, la quasi-totalité des atomes donneurs sont ionisés. 

Si N d est la concentration des atomes donneurs, ceux-ci vont libérer n = N d électrons 
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Les concentrations en électrons libres (n) et en trous libres (p) sont liées par la loi d’action 
de masse : 

n.p = nf 1. 3 

Les électrons sont les porteurs majoritaires et les trous les porteurs minoritaires. 

Le niveau indicateur de Fermi E Fji se déplace donc du milieu de la bande interdite (E F .) vers 
la bande de conduction de telle manière que : 

/N d \ 

A E n = KT. In I . 4 

Avec : A E n = E Fn - E F . 

1.4.2 Semi-conducteurs du type p 

On obtient un S.C. dopé p en injectant dans le silicium des atomes de la 3° colonne (bore, 
indium) qui possèdent trois électrons périphériques (figure 1.3). 




O 

& 


Bande de conduction 



Ev 


Bande de valence 



niveau de Fcnni EFp 


\ 


A£„ = ÀTln(— ) 
” w, 


P S K i 


Fig.I.3 : Libération d’un trou par l’atome de bore et schéma des bandes. 

A la température ordinaire, la quasi-totalité des atomes accepteurs sont ionisés. Si N a est la 
concentration par cm' 3 des atomes accepteurs, ceux-ci vont libérer :p = N a trous libres. 

Les concentrations en électrons libres (n) et en trous libres (p) sont liées par la loi d’action 
de masse : 

n.p = nf 
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Les trous sont les porteurs majoritaires et les électrons les porteurs minoritaires. 

Le niveau indicateur de Fermi E Fp se déplace du niveau intrinsèque E F . vers la bande de 

valence de telle manière que : 

A E p = KT. In 1. 6 

Avec : A E p = E F . - E Fp 


1.5 Propriétés du matériau CIGS 

Le CIGS est le matériau clé dans les cellules CIGS dont il convient d’étudier les propriétés 
structurales et optoélectroniques. 

1.5.1 Structure de CIGS 


Le CIGS est un acronyme de séléniure de cuivre-indium gallium (CuIn 1 _ x Ga x Se 2 ) est 
un semi-conducteur quaternaire en solution solide ou x représente la portion du gallium. La 
concentration d'indium et de gallium peut varier du séléniure de cuivre et d'indium pur (x=0, 
CIS) à du séléniure de cuivre et de gallium pur (x=l, CGS). Ainsi, son gap varie 
continuellement avec x, de 1.0 eV pour CIS au environ de 1.7 eV pour le CGS. 
Cristallographiquement, le CIGS est un tétraèdre appartenant à la famille des I-III- VI 2 avec une 
structure de chalcopyrite (Fig.1.4). 



• II- Cu 


# m-In/Ga 
O VI Se/S 



Cu 

• lu 

• Se 

• Ga 


Fig.1.4 : Structure cristalline de la chalcopyrite [1,2]. 
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Les paramètres de maille sont contrôlés par la composition du ternaire (CuInSe 2 , CuGaSe 2 , 
CuInS 2 , CuGaS 2 ...), ils suivent la loi de Vegard [3-4], Pour les chalcopyrites intégrant les 
atomes les plus petits (Ga au lieu de In, S au lieu de Se), la maille cristalline se contracte, donc 
le paramètre de maille a diminué (Fig.1.5), et inversement pour les atomes les plus gros. 

Les variations de composition ont également une influence sur les valeurs de bande interdite, 
les énergies des bandes interdites varient entre 1,0 eV pour CuInSe 2 et 2,4 eV pour CuGaS 2 
[5]. D’une façon tout à fait qualitative, on peut dire que pour des atomes plus lourds avec une 
polarisabilité plus grande (In par rapport au Ga ou Se par rapport au S), les bandes de 
valence/conduction s’élargissent, et la valeur de la bande interdite de la chalcopyrite diminue. 

La figure 1.5 représente la valeur de la bande interdite en fonction de la valeur du paramètre 
de maille a, pour différentes chalcopyrites. 

Ce diagramme montre que la largeur de la bande interdite varie entre 1 et 1 .68 eV pour les 
deux composés CuInSe 2 et CuGaSe 2 respectivement en variant la composition en Ga dans 
l’alliage CuQn^x, Ga x )Se 2 . 
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Fig.1.5 : Largeurs de bande interdite et paramètre de maille a de certaines chalcopyrites 
d’après [5], 
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1.5.2 Ajustement de la bande d’énergie interdite 

Comme indiqué précédemment, le CIGS est un alliage entre le CIS et le CGS. La 
représentation de CIGS est : 

CuÇln^x, Ga x )Se 2 1. 7 

Où x représente la fraction molaire de Ga qui est présente dans l'alliage et donnée par : 

Gu 

x = — — — 1. 8 

(In + Ga) 

Une valeur x=0 représente le CIS pur avec une bande interdite de 1,07 eV. Une valeur x=l 
représente le CGS pur avec une bande interdite de 1,76 eV [6]. Pour x variant de zéro à un, le 
fabriquant peut ajuster la bande interdite du CIGS à la valeur désirée. L'ajustement de la bande 
interdite selon le contenu de Ga ne suit pas une loi linéaire. L'équation reliant de bande interdite 
à la proportion x de gallium pour les CIGS intrinsèques à 300 K est comme suit [7] : 

E g = 1.011 + 0.664x - 0.249x(l - x) 1. 9 


A partir de (Equ.1.9), la relation entre x et la bande interdite est légèrement parabolique. 
L'énergie de la bande interdite est liée à la longueur d'onde du photon par [6] : 


E g (eV) = 


1.24 

A(/rm) 


1.10 


Où E g est l'énergie de la bande interdite en électrons volts et X est la longueur d'onde de la 
lumière en micromètres. D'après (Equ.1.10), la bande interdite est inversement proportionnelle 
à la longueur d’onde. Pour le matériau à large bande interdite, il ne peut absorber efficacement 
que la lumière de longueurs d'onde courtes. 

1.5.3 Propriétés optiques 

Le coefficient d’absorption de CuInSe 2 est très élevé et de l’ordre de 10 5 cm' 1 [8] dans le 
domaine du visible et du proche infrarouge (Fig.1.6). 

L'équation de conversion des longueurs d'ondes du spectre solaire en énergie est la suivante : 
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Avec : 

E : l'énergie d'irradiation en eV. 

h : la constante de Planck, elle est égale à 6.62 xlO' 34 J.s. 
c : la vitesse de la lumière dans le vide, elle est égale à 3xl0 8 m s" 1 . 

X : la longueur d'onde incident en m. 

La valeur de l’énergie du gap du CuInSe 2 permet d'absorber les photons ayant des longueurs 
d'onde inférieures à 1192 nm, ce qui indique que la gamme d'absorption de ce ternaire sont 
située dans la région optimale du spectre solaire. 



Fig.1.6 : Coefficient d’absorption des différents matériaux [9]. 


Rappelons que le CuInSe 2 est un semiconducteur à gap direct, la relation entre son gap 
d’énergie E g et son coefficient d'absorption a est la suivante [10] : 


a = 


A(E — Eg)^ 


14 


Y 


E 
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Avec A, une constante de proportionnalité qui dépend des densités des états liés à 
l'absorption de photon, E l'énergie d'irradiation. 

La figure 1.7 représente le coefficient d’absorption d’une couche mince en CuInGaSe 2 avec 
x = Ga/(In + Ga) = 0 et 0.3 en fonction de longueur d’onde. 



Fig.1.7 : Coefficient d'absorption de CuInSe 2 et CuInGaSe 2 en fonction de longueur 
d’onde [10]. 


Le gap d’énergie de ces semiconducteurs est également une fonction de la température 
(Fig.1.8), il peut être décrit approximativement par la relation proposée par Varshni [1 1] : 


Eg(T) = Eg( 0) - 



1.13 


Où Eg (0) est le gap à 0 K, a et /? sont des constantes spécifiques à chaque semi-conducteur. 
Dans le cas du CuInSe 2 , ils sont respectivement égaux à 3.6 xlO 4 eV. K' 1 et 350 K [12], 
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Fig.1.8 : Variation de l’énergie du gap de CuInSe 2 en fonction de la température [12]. 


Notons aussi la relation suivante entre le coefficient d’extinction et l’indice de refraction : 

N = n + ik 1. 14 

Où N représente l’indice optique complexe, k est le coefficient d’extinction du matériau 
traduisant l’absorption du rayonnement par le matériau et n est l’indice de réfraction. C’est le 
rapport entre la vitesse de propagation de la lumière dans le vide et la vitesse dans le milieu 
étudié. 

Le coefficient d’extinction k est lie au coefficient d’absorption a, pour chaque longueur 
d’onde k, par la relation : 


a = 


Ank 

~T 


1.15 
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La figure 1.9 représente l’évolution de l’indice de réfraction n et le coefficient d’extinction 
k du CIGS en fonction de la longueur d’onde. 



Fig.1.9 : L’évolution de l’indice de réfraction n et le coefficient d’extinction k du CIGS en 
fonction de la longueur d’onde [13]. 

1.6 Les techniques d’élaboration du matériau CIGS 

Les techniques actuelles de préparation de CIGS se divisent en trois grandes catégories : la 
coévaporation, le procédé de sélénisation/sulfuration de précurseur métallique et les procédés 
de fabrication atmosphérique. Le contact arrière (c’est-à-dire la couche de molybdène) est 
déposé par pulvérisation cathodique. La sélénisation (l’ électrodéposition ou la pulvérisation 
suivie d’une exposition à un gaz contenant du sélénium), la coévaporation ou la pulvérisation 
réactive sont ensuite utilisées pour déposer la couche absorbante. La strate tampon 
naturellement dopée n est adjointe via un bain chimique. Enfin, la couche fenêtre (conducteur 
supérieur transparent) est toujours déposée par pulvérisation. La meilleure cellule CIGS a été 
fabriquée par coévaporation [14], 
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1.6.1 Evaporation à partir des sources élémentaires 

La coévaporation est une méthode de dépôt physique par phase vapeur (ou PVD : Physical 
Vapor Déposition). Les matériaux purs à déposer sont placés sous forme solide dans des 
creusets portés à haute température. La température des creusets conditionne les flux de vapeur 
de chaque élément et donc la composition du dépôt. Un vide est nécessaire pour permettre aux 
flux de vapeur d’atteindre le substrat et pour éviter les impuretés. Le substrat est porté à une 
température permettant la diffusion des éléments et le contrôle de la croissance cristalline du 
matériau au fur et à mesure de l’arrivée des vapeurs réactives [15]. Un schéma de principe de 
cette technique est présenté dans la Figure 1. 10 suivante : 



Fig.1.10 : Principe de la technique de la coévaporation thermique pour la couche CIGS [14]. 

Cette méthode conduit à une marge de manœuvre aisée sur le rapport de Ga/(Ga+In). Ce 
dernier permet d’obtenir un diagramme de bande mieux adapté au spectre de rayonnement 
incident. Très récemment, grâce à cette technique, des rendements élevés ont été réalisés par 
ZSW (20.8%) et NREL (20.3%) [16]. Le procédé de coévaporation à partir de sources 
élémentaires, optimisé pour obtenir des cellules de haut rendement, se compose typiquement 
de trois étapes [17], présentées dans la figure 1.11. 
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Fig.I.ll : Flux de métaux et profil de température de substrat utilisés dans un procédé en 
trois étapes de synthèse du CIGS par coévaporation [18]. 

Une première étape implique le dépôt de (In,Ga) 2 Se 3 sur un substrat à 300-350°C. Une 
deuxième étape de dépôt de Cu et Se permet d’obtenir une couche riche en cuivre et d’initier 
la croissance de gros grains, alors que le procédé se termine par un dépôt riche en indium et 
gallium pour assurer les bonnes propriétés électroniques de la couche. Les dernières étapes sont 
effectuées avec une température de substrat de 550 à 600°C. Le refroidissement se fait sous 
vapeur de Se [19]. L’avantage de cette technique est de pouvoir ajuster facilement la 
composition, et d’obtenir des films de haute cristallinité avec peu d’impuretés. Mais le 
changement d’échelle pour obtenir de grands panneaux reste difficile, et la technique est 
coûteuse en énergie à cause du vide poussé (de l’ordre de 10" 8 mbar). De plus, elle ne permet 
qu’une utilisation partielle des matériaux utilisés dans les creusets d’évaporation puisque le 
matériau se dépose aussi sur les parois de l’évaporateur [15]. 

1.7 Conclusion 

Dans ce chapitre nous avons présenté quelques généralités de base sur la physique des semi- 
conducteurs et quelques propriétés du matériau CIGS utilisé dans la conception de la cellule 
solaire à étudier dans ce travail. Ces propriétés concernent les propriétés structurales, 
électroniques, électriques et optiques. Nous avons aussi décrit le procédé d’élaboration du 
CIGS. 
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II. 1 Introduction 

La conversion photovoltaïque est la transformation directe d’une énergie électromagnétique 
(rayonnement) en énergie électrique au moyen d’une cellule solaire, l’effet photovoltaïque a 
été utilisé au début des années cinquante dans les applications pour l’espace (alimentation des 
satellites) et ses premières utilisations au niveau de la terre est dans les années soixante-dix. 

L’élément de base dans la conversion du rayonnement de la lumière est la cellule 
photovoltaïque. Plusieurs cellules sont associées dans un module. Un ensemble de modules 
sont groupés en panneaux eux même associés pour former le champ photovoltaïque. 

Les cellules solaires font actuellement l’objet de multiples recherches dans le but de réaliser 
le meilleur rapport entre le rendement énergétique et le prix de revient. 

Dans ce chapitre, nous débuterons par donner les caractéristiques du rayonnement solaire et 
la description de la conversion photovoltaïque tout en étudiant les phénomènes physiques ayant 
lieu tels la génération des photoporteurs et leur recombinaison dans un semiconducteur. 

Ensuite, nous décrivons le principe de fonctionnement et les caractéristiques des cellules 
solaires en général et les cellules solaires en CIGS en particulier. 


II.2 Rayonnement solaire 

La distance de la terre au soleil est d'environ 150 millions de kilomètres et la vitesse de la 
lumière est d’un peu plus de 3xl0 8 m/s, les rayons du soleil mettent donc environ 8 min à nous 
parvenir [1], 

Le rayonnement reçu par la terre en incidence normale avoisine les 1400 W/m 2 avant 
d’aborder la couche atmosphérique [1]. 

Pour tenir compte de l’épaisseur d’atmosphère traversée par le rayonnement solaire incident 
(Ligure II. 1), on introduit un coefficient appelé nombre de masse (AM) défini par : 


AM = 


1 

cos(9) 


II. 1 


Où 0 est l’angle que fait le soleil avec son zénith. Par définition, AMI correspond aux 
Conditions hors atmosphère. 

Quand le soleil fait un angle de 48° par rapport au zénith, la lumière incidente est dite 
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Le nombre de masse d'air caractérise la puissance transportée par le rayonnement solaire 
(83.3 mW/cm 2 pour AM 1.5) et de plus, sert de standard pour quantifier les performances de 
nouveaux dispositifs photo voltaïques. 

Les conditions standard de caractérisation sont définies par les normes IEC-60904 de 
l'International Electrotechnical Commission (IEC) selon une distribution spectrale AM1.5G 
global (somme des rayonnements directs et diffus), d’intensité 100 mW/cm 2 et pour une 
température de cellule de 25 °C [1]. 

Au sol, le rayonnement solaire a au moins deux composantes : une composante directe et 
une composant diffusée (rayonnement incident diffusé ou réfléchi par un obstacle : nuages, sol) 
formant le rayonnement global. 

La figure II. 1 montre le spectre d’émission solaire sous AM 1.5G global normalisé à 
100 mW/cm 2 [1], 


le soleil au zénith 





Fig.II.l : Radiation spectral pour différent air masses. 


Le spectre solaire c'est sa décomposition en longueurs d'onde ou "couleurs". La lumière 
solaire est en effet composée de toutes sortes de rayonnements de couleurs différentes, 
caractérisées par leur gamme de longueur d'onde. Les photons, grains de lumière qui composent 
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ce rayonnement électromagnétique, elles sont porteuses d'une énergie qui est reliée à leur 
longueur d'onde par la relation : 

hc 1.24 

E = hv = — = — — II. 2 

À À 

Où v : fréquence de la radiation, X : longueur d’onde de la radiation en pm, c : vitesse de la 
lumière, E : énergie du photon en eV, et h : constante de Planck. 

La figure IL 2 représente le spectre solaire et la réponse spectrale des matériaux en CIGS. 



Fig.II.2 : Spectre solaire et la réponse des matériaux CIGS [2], 

II.3 Énergie solaire photovoltaïque 

L'énergie solaire photo voltaïque est une énergie électrique produite à partir du rayonnement 
solaire. L'énergie produite par une centrale solaire photovoltaïque est dite renouvelable, car sa 
source est considérée comme inépuisable à l'échelle de temps humaine. Cette énergie n'est 
cependant pas totalement renouvelable, car la fabrication des panneaux photovoltaïques, leur 
installation et leur exploitation consomment de l'énergie issue en grande partie de sources non 
renouvelables ; mais un système photovoltaïque produit entre 20 et 40 fois plus d'énergie que 
ce qui a été utilisé pour le fabriquer [3]. 
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La cellule photovoltaïque est le composant électronique de base. Elle utilise l'effet 
photoélectrique pour convertir les ondes électromagnétiques émises par le soleil en électricité. 
Plusieurs cellules reliées entre elles forment un module solaire photovoltaïque, plusieurs 
modules regroupés forment une installation solaire. L'électricité est soit consommée ou stockée 
sur place, soit transportée par le réseau de distribution et le réseau de transport. 

En novembre 2015 a été lancée l'alliance internationale pour l'énergie solaire (ou ISA pour 
« International solar alliance »), projet porté par Narendra Modi, visant à réunir les états riches 
en ressources solaires pour coordonner les politiques de développement du solaire thermique 
et photovoltaïque (formation, développement de standards communs, partage d'expériences, 
coentreprises...). La cérémonie de lancement, organisée par l'Inde et la France, a eu lieu à Paris, 
avec le soutien du secteur privé [3]. 

En 2016 la consommation mondiale d'électricité s'élève à 23 816 TWh/an, soit l'équivalent 
de la production d'une surface photovoltaïque de 100 000 km 2 (équivalent de la superficie de 
l'Islande) [3], 

II.4 Génération de porteurs en excès par la lumière 

Une radiation lumineuse rencontrant un semi-conducteur est absorbée suivant la loi de 
Lambert-Bouguer [4] : 

I(x) = 7(1 — R) exp(—ax ) = 7 0 exp(—ax ) II. 3 

Où I : le flux de photons (cm'V 1 ), x : profondeur d’absorption du faisceau dans le matériau 
à partir de la surface du semiconducteur, R : coefficient de réflexion, représente la part de 
l’énergie lumineuse incidente /, réfléchie à la surface du matériau, a : coefficient d’absorption, 
traduit la probabilité d’absorption d’un photon par unité de longueur. 

Quand les photons incidents ont une énergie suffisante pour provoquer le transfert des 
électrons de la bande de valence vers la bande de conduction, des paires électron-trou peuvent 
se créer. Ce phénomène appelé mode d’absorption intrinsèque ou fondamentale prédomine 
dans les cellules solaires [4], 

Le calcul du taux de génération de paires électron-trou, G (x) (cm'V 1 ) s'effectue de la façon 
suivante. Appelant 7 0 , le flux de photons (cm'V 1 ) incidents sur la face éclairée du matériau et 
a (cm.- 1 ) le coefficient d'absorption de la lumière par le matériau. A la profondeur x, ce taux 
de génération vaut [4] : 
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G(x) = a. (1 — R)I 0 . exp(— ax) II. 4 

II.5 Mécanismes de recombinaison 

Nous avons vu que dans une cellule solaire, des paires électron-trou sont créées par 
absorption de la lumière. Les recombinaisons de ces paires électron-trou font partie des 
phénomènes fondamentalement limitants pour la performance des dispositifs photovoltaïques. 

Un mécanisme de recombinaison donné est caractérisé par un taux de recombinaison R, 
représentant le nombre de recombinaisons par unité de temps et par unité de volume de 
matériau semiconducteur. Il existe trois principaux mécanismes : l’effet Auger, la 
recombinaison radiative et la recombinaison assistée par un piège dans la bande interdite. Cette 
dernière est aussi appelée recombinaison Shockley-Read-Hall (SRH). Ce dernier mécanisme 
est le dominant dans les cellules basées sur des matériaux quaternaires et polycristallins comme 
le CIGS. La figure II.3 illustre ces différents mécanismes au sein d’un semiconducteur [5,6]. 



Fig.II.3 : Principaux mécanismes de recombinaison dans un semiconducteur : 
Shockley-Read-Hall (SRH), radiative et Auger [6,7], 
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Les recombinaisons SRH se produisent lorsque des niveaux de pièges présents dans la bande 
interdite où se manifeste la capture des porteurs. Les niveaux susceptibles de participer aux 
recombinaisons SRH peuvent provenir de défauts structurels ou de la présence d’impuretés 
dans le CIGS. Le taux de recombinaison par ce mécanisme pour un piège d’énergie E p vaut 
[6,7]: 

pn — nj 

RsRH = “ Et-E p \ , ~ E p - E t \ 11 5 

*sRH,n (P + n t exp y J + t SRH)P + n^exp ^ J 


Avec t srh n et t srh , p les temps de vie de chaque type de porteurs et E t le niveau de Fermi 
dans le semiconducteur intrinsèque. Cependant, il est possible de simplifier cette expression. 
En effet par définition dans un semiconducteur dopé, les porteurs majoritaires sont en large 
excès par rapport aux porteurs minoritaires. Puis qu’il est nécessaire de disposer d’un porteur 
de chaque type pour qu’une recombinaison ait lieu, c’est donc la densité de porteurs 
minoritaires qui détermine principalement le taux de recombinaison. Ainsi, dans des conditions 
normales de fonctionnement d’une cellule, pour un semiconducteur de type p comme le CIGS 
et en considérant que les niveaux de piège sont proches du milieu de la bande interdite 

(j E p ~ ^ le taux de recombinaisons SRH peut s’écrire [6-7] : 


R 


SRH ~ 


n — n 0 

T SRH,n 


II. 6 


Avec n 0 la densité d’électrons à l’équilibre. Le temps de vie des électrons peut s’exprimer 
de la façon suivante : 

1 

T SRH,n ~ —— TT H- ? 

^■nVth^p 

Avec <r n la section efficace de capture des électrons par le piège considéré, v th la vitesse 
thermique des électrons et N p la densité de pièges d’énergie E p . Ainsi, le temps de vie des 
porteurs sera réduit si la densité de pièges augmente. Il est donc évident qu’un matériau 
présentant plus de défauts structurels ou d’impuretés présentera un taux de recombinaisons 
SRH grand. Ce mécanisme est dominant dans les cellules basées sur des matériaux quaternaires 
et poly cristallins comme le CIGS. Les expressions précédentes étant valables pour un piège 
avec une énergie donnée, le taux total de recombinaisons SRH vaut [6] : 
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R 


SRH 


-I 


R 


SRH,i 


II. 8 


Chaque i étant un piège avec un niveau d’énergie E p i dans la bande interdite. 

Les recombinaisons radiatives se font par transition directe d’un électron de la bande de 
conduction vers la bande de valence. L’énergie de la transition est libérée sous forme d’un 
photon. Dans le CIGS de type p, en tenant compte des mêmes simplifications que pour les 
recombinaisons SRH, le taux de recombinaisons radiatives s’écrit [6] : 


n — Un 

D ~ 

lx rad ~ 

T Rad,n 


II. 9 


Le temps de vie des électrons T Rad n vis-à-vis des recombinaisons radiatives vaut[6] : 


1 


T-RacL.n ~ 


P oB 


II. 10 


Avec B le coefficient de recombinaison radiative et p 0 la densité de trous à l’équilibre 
thermodynamique . 

On remarque donc que pour les recombinaisons radiatives, le temps de vie des porteurs 
diminuera quand le dopage p du CIGS augmente. Le photon émis lors d’une recombinaison de 
ce type possède une énergie proche du gap du matériau. Il est donc possible que celui-ci soit 
réabsorbé et reforme ainsi une paire électron-trou. Cette récupération n’est pas possible dans 
le cas d’une recombinaison SRH [6], 

Les recombinaisons Auger sont elles aussi des transferts directs bande à bande d’un porteur. 
Dans ce cas il peut s’agir d’un électron ou d’un trou. Cependant, au lieu d’être émise sous 
forme d’un photon, l’énergie est transférée à un autre porteur du même type. Ce dernier se 
thématisera ensuite par émission de phonons [6] (voir figure II.3). 

Pour le CIGS, dans le cas simplifié, le taux de recombinaisons Auger s’écrit [6] : 

n — Tin 

P ~ 

lx Auger ~ 

T Auger, n 

Le temps de vie des électrons T Auger n vis-à-vis des recombinaisons Auger vaut [6] : 

1 

T Auger,n ~ ^2^ 
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Où C p est le coefficient de recombinaison Auger. On voit ici que le temps de vie des 
porteurs diminue suivant le carré de la densité de trous et donc du dopage. Cela montre que 
les bénéfices d’un dopage plus élevé sont limités par les recombinaisons Auger. 

Des niveaux d’énergie d’états spécifiques (états de surface), provenant de phénomènes 
intrinsèques (des liaisons non satisfaites des atomes de surface, disparition de la périodicité du 
réseau, adsorption d’atomes étrangers), peuvent se situer dans la bande interdite. Certains de 
ces états jouent également le rôle de centres de recombinaison. La surface d'un semiconducteur 
est un lieu où les recombinaisons sont nombreuses. Le taux de recombinaison Shockley-Read- 
Hall en surface est donné par l’expression suivante [8] : 


U s 


n sPs - nj 

Oh±rh) + (Ps + Pô 


11.13 


Où S p et S n (cm/s) sont respectivement les vitesses de recombinaison en surface pour les 
électrons et les trous, n s et p s sont les concentrations des porteurs à la surface. et p x 
représentent les termes d'émission et sont donnés par [8]: 


ni 


riiexp 


E jzEl± 

KT 


Px = r^exp 


-(gt-g/ f ) 

KT 


II. 14 


Tels que E t représente l'énergie du niveau piège proche du milieu de la bande interdite, Ef. 
le niveau de Fermi intrinsèque et ni la concentration intrinsèque. 

Pour mieux comprendre le principe de fonctionnement d’une cellule photo voltaïque, il est 
important de décrire les phénomènes qui régissent l’interaction entre un semi-conducteur et le 
rayonnement lumineux. 

II.6 Interaction rayonnement/semi-conducteur 

Les trois processus qui peuvent se produire lors de l’interaction rayonnement/semi- 
conducteur sont [9] : 

4- L’absorption fondamentale qui représente le passage d’un électron d’un état occupé 
de la bande de valence à un état vide de la bande de conduction après l’absorption 
d’un photon, 

4- L’émission spontanée qui représente le passage d’un électron de la bande de 
conduction à la bande de valence avec émission d’un photon, 
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■4- L’émission stimulée qui représente le passage d’un électron de la bande de conduction 
à la bande de valence suite à l’absorption d’un photon avec l’émission d’un deuxième 
photon. 

Ces différents processus sont fonction du type de semi-conducteur. Nous devons distinguer 
les semi-conducteurs à bande interdite directe et ceux à bande interdite indirecte (Fig.II.4). 


a) 


hu A 

1 Absorption 



b) 



D'un pliaiii 


Fig.II.4 : Transitions bande à bande dans un semiconducteur (a) à bande interdite directe 

et (b) à bande interdite indirecte [10], 

Un semi-conducteur est à bande interdite directe quand le minimum de la bande de 
conduction et le maximum de la bande de valence sont verticales (voir figure II.4(a)). Les 
transitions bande à bande (de la bande de conduction vers la bande de valence) sont radiatives 
(émission spontanée) et sont favorisées avec ce type de semi -conducteur utilisé en 
optoélectronique (matériau III- V en général comme le GaAs) [9]. 

Un semi-conducteur est à bande interdite indirecte quand le minimum de la bande de 
conduction est décalé par rapport au maximum de la bande de valence. 

Les transitions bande à bande sont non-radiatives (voir figure II.4(b)). Ceci est le cas même 
si avec un photon de plus grande énergie, un électron peut être directement excité vers un 
minimum relatif de la bande de conduction et se thermalisme dans le minimum absolu de la 
bande de conduction. La transition d’un électron dans un matériau à bande interdite indirecte 
est assistée par un phonon absorbé ou émis par l’électron afin que son vecteur d’onde 
corresponde au maximum de la bande de valence pour absorber un photon [9]. 
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II.7 Principe de fonctionnement d’une cellule photovoltaïque 

Le fonctionnement de toute cellule photovoltaïque est basé sur l’existence d’un champ 
interne afin de séparer les paires électron- trou générées par l’absorption de la lumière. 

Cette barrière de potentiel peut être créée par [9] : 

4- Une jonction p-n avec un même semi-conducteur que l’on appelle homojonction 
■4- Une jonction p-n avec deux semi-conducteurs différents que l’on appelle 

héîéroj onction 

4- Un contact métal/semi-conducteur de type Schottky. 

Le principe de fonctionnement d’une cellule à jonction p-n à homojonction est décrit dans 
la figure IL 5. 


N 

1 


P 

JL_ 





Fig.II.5 : Structure et diagramme de bande d’une cellule solaire. 


Les photons incidents créent des paires électron-trou dans les régions N et P et dans la zone 
de charge d’espace [9-11] : 

4- Si une paire est générée dans une zone électriquement neutre (P ou N), les porteurs 
diffusent. Si les porteurs minoritaires atteignent la zone de charge d’espace alors ils 
sont propulsés par le champ électrique dans la zone où ils deviennent porteurs 
majoritaires. Ces porteurs contribuent ainsi au courant par leur diffusion et il se crée 
un courant de diffusion. 

4- Si une paire électron-trou est générée dans la zone de charge d’espace, l’électron et le 
trou sont séparés par le champ électrique et chacun est propulsé dans la région où il 
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est porteur majoritaire (région N pour l’électron et région P pour le trou). Ces porteurs 
donnent naissance à un courant de génération. 

Ces deux contributions donnent naissance à un photocourant qui contribue au courant 
inverse de la diode formée par la jonction. Le courant inverse d’une jonction p-n est fonction 
des densités de porteurs minoritaires dans les régions neutres de la diode et de la génération de 
paires électron-trou dans la zone de charge d’espace. L’éclairement augmente le courant 
inverse par la création de porteurs minoritaires dans les zones neutres et la génération de paires 
électron-trou dans la zone de charge d’espace [9]. 

II.8 Caractéristique I (V) et schéma équivalent d’une cellule solaire 

Une cellule solaire reliée à une résistance électrique et soumise à l'éclairement solaire débite 
un courant qui va se répartir entre la résistance et sa structure interne de diode. Sa 
caractéristique I (V) correspond à la soustraction du photocourant et du courant de la diode à 
l’obscurité par [12] : 


m = un - v = h ( ex P (|£) - 1) - i n ». 15 

Avec I ph le photocourant, I obs le courant d'obscurité, I s le courant de saturation de la diode, 
q la charge élémentaire, K la constante de Boltzmann et T la température. 

La caractéristique d'une cellule à l'obscurité est identique à celle d'une diode. Nous avons 
représenté sur la figure IL 6 les deux caractéristiques courant-tension d’une cellule solaire à 
l’obscurité en pointillés, et sous illumination en trait plein. 

On remarque que la courbe sous illumination est simplement décalée par rapport à la 
première d’une valeur I cc , qui traduit la génération constante du courant par la lumière. Cette 
valeur est appelée courant de court-circuit, puisque c’est le courant que génère la cellule sous 
lumière à tension nulle. La valeur V co , à l'inverse, est la tension de circuit ouvert. Il représente 
en fait l’égalité du courant de génération par la lumière et du courant à l'obscurité [13]. 
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I 



Fig.II.6 : Caractéristiques I(V) à l'obscurité et sous éclairement d’une cellule photo voltaïque. 

Son schéma équivalent est représenté par une diode idéale montée en parallèle avec une 
source de courant (figure II.7). La résistance série, R s , cette résistance est liée à l’impédance 
des électrodes et de la base, il en résulte que la tension V aux bornes de la cellule est différente 
de la tension aux bornes de la jonction. La résistance shunt, R sh , qui correspond aux pertes 
dans la surface et pertes dues aux défauts dans le matériau, il en résulte qu’une partie du courant 
I ph sera dérivée par cette résistance et ne pourra être délivrée à la charge [14]. 


R s 



Fig.II.7 : Schéma électrique réel d’une cellule photovoltaïque. 


Le courant / peut s’écrire : 
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Selon le schéma équivalent d’une cellule solaire dans la figure II.7, on a : 



11.17 



II. 18 



II.9 Paramètres d'une cellule solaire 

Les grandeurs caractéristiques principales des cellules solaires sont : 

■J- Le courant de court-circuit I cc 
4- La tension en circuit ouvert V co 
4- Le facteur de forme FF 
4- Le rendement de conversion d’énergie rj 
4- La réponse spectrale RS 

11.9.1 Courant de court-circuit, I cc 

Le courant de court-circuit est le courant débuté par la cellule sous éclairement en court- 
circuitant la sortie. C’est à dire que : I cc = I(V = 0). 

Pour une cellule solaire idéale le courant de court-circuit est égal au courant de photo 
génération I ph [14], 

11.9.2 Tension de circuit ouvert, V co 

La tension de circuit ouvert est la différence de potentiel mesurée aux bornes de la cellule 
lorsque la sortie est en circuit ouvert. C’est à dire : V co = V (/ = 0). 

Pour une cellule solaire idéale la tension de circuit ouvert est donnée par [14] : 
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II.9.3 Facteur de forme, FF 

La puissance fournie au circuit extérieur par une cellule photovoltaïque sous éclairement 
dépend de la résistance de charge (résistance externe placée aux bornes de la cellule). Cette 
puissance est maximale (P max ) P our un point de fonctionnement P ma xOmax« V max ) de la courbe 
courant-tension (Figure II.6). Le facteur de forme est défini par la relation suivante [13] : 

P VI 

r max v max- l max ¥¥ 

FF= n — n~ 1120 

v CO ml CC v CO' 1 CC 


II.9.4 Rendement de conversion, // 

Le rendement rj, des cellules solaires désigne le rendement de conversion en puissance. Il 
est défini comme étant le rapport entre la puissance maximale délivrée par la cellule et la 
puissance lumineuse incidente P in [ 13]. 


P VI 

r max rp v ccr L cc 

77 ~ — = FF. 

P. P. 

1 in 1 in 


11.21 


Où P in est la puissance incidente de la lumière est égale à la puissance solaire P so iaire 
C Psolaire = 100 mW/cm 2 ). 

Ce rendement peut être amélioré en augmentant le facteur de forme, le courant de court- 
circuit et la tension en circuit ouvert. 


A température et éclairage constants, le rendement d’une cellule solaire dépend de la charge 
dans le circuit électrique. En circuit ouvert ( R c = oo, 7 = 0, V = V C0 ) ou en court-circuit 
(7? c = 0, I = l cc , V = 0), aucune énergie n’est transmise à l’extérieur. Entre ces deux 
extrêmes, il existe une valeur optimale R opt de la résistance de charge R c pour laquelle la 
puissance fournie par la cellule solaire à la résistance de charge est maximale notée P max [13] 
(figure II. 8). 
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Fig.II.8 : Caractéristiques I(V) et P(V) d'une cellule solaire. 

II.9.5 Réponse spectrale d’une cellule photovoltaïque 

La réponse spectrale RS d'une cellule photo voltaïque est le rapport entre le courant de court- 
circuit généré par la cellule et la puissance lumineuse incidente, en fonction des différentes 
longueurs d'onde formant le rayonnement incident. 

La réponse spectrale RS est donnée par la relation suivante [15] : 

I C M l) 

RSW= R(é (A/W) 1122 

Le rendement quantique externe EQE de la cellule est le rapport du nombre de porteurs 
générés sur le nombre de photons incidents pour chaque longueur d'onde, il est relié à la réponse 
spectrale par [15] : 


EQE(X) = RS(X)— 11.23 

Aq 

Le rendement quantique interne 1QE de la cellule est le rapport du nombre de porteurs 
générés sur le nombre de photons pénétrant dans le matériau, il est relié au rendement quantique 
externe EQE par l'équation [15] : 


IQE(A ) = 


EQE (A) 

1 -RW 


11.24 
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Le coefficient de réflexion s’écrit 


RW = ( 


n 2 (A)-n 1 (A) 
n 2 U)+n!(/ 1) 



matériaux en contact sont 


d’indices de réfraction n x et n 2 . 


II. 10 Cellule solaire en couche-mince 

Par principe, une couche mince est une fine couche d’un matériau déposé sur un autre 
matériau, appelé "substrat" dont l’une des dimensions qu’on appelle l’épaisseur a été fortement 
réduite de telle sorte qu’elle varie de quelques "nm" à quelques "pm" (typiquement ce sont des 
couches de 10 ... 100 nanomètres d'épaisseur). 

Dans le domaine des couches minces, il existe trois principales filières : 

’L La filière silicium amorphe et microcristallin noté TFSi (Thin-Film Silicon). 

■4- La filière CdTe (Tellurure de Cadmium). 

=J- La filière Cu(In, Ga)Se 2 ( Cuivre-Indium/Gallium-Sélénium ) noté CIGS. 

II. 10.1 Filière silicium en couches minces 

Depuis les années 1976, des recherches intensives ont été entreprises pour utiliser du 
silicium non cristallisé, c’est-à-dire à l’état amorphe, dont les liaisons pendantes sont saturées 
par hydrogénation (a-Si :H). Ce matériau présente un fort coefficient d’absorption qui autorise 
de très faibles épaisseurs de l’ordre du micron et une faible consommation énergétique durant 
le cycle de production et l’aptitude à être déposé en grandes surfaces. Par contre, ses deux 
points faibles sont le rendement de conversion et la dégradation sous lumière (instabilité dite 
de Staebler-Wronski). 
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La figure II.9 présente la plus simple structure d'une cellule en silicium amorphe, elle se 
compose d'une zone dopée bore, intrinsèque et dopée phosphore (p-i-n). 


Silicium amorphe 
Lumière 


Verre 

Conducteur 

transparent 

Couche p 
Couche / 
Couche n 
Métal 



+ 


Fig.II.9 : Cellule solaire au silicium amorphe [16]. 

On peut trouver dans la figure II. 10, les caractéristiques électriques réelles et simulées d’une 
cellule en silicium amorphe développée en laboratoire [17]. Cette cellule montre un rendement 
record de 8.5% pour une simple jonction amorphe avec un facteur de forme de 0,719. On peut 
remarquer l’existante d’une différence importante des performances entre les différentes 
technologies silicium. La cellule amorphe produit 22.5mA/cm 2 contre 38.1mA/cm 2 pour une 
cellule solaire en silicium multi-cristalline avec une tension de circuit-ouvert de 0.524V contre 
0.654V de la cellule solaire en silicium multi-cristallin [18]. 



Fig.II.10 : Caractéristique courant-tension d’une cellule solaire en silicium amorphe, (a) 
Courbe réelle, (b) Courbe simulée [18]. 
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II.10.2 Filière tellurure de cadmium (CdTe) 

Les cellules solaires au tellurure de cadmium (CdTe) présentent les avantages liés à la 
technologie en couches minces, elles sont légères et robustes. Le tellurure de cadmium est 
employé en tant que semi-conducteur dopé p. Un autre matériau est donc requis pour réaliser 
la jonction p-n, à savoir du sulfure de cadmium dopé n (CdS) (Fig.II.ll). 


Verre 


Coûta et a vaut 
CdS 

CdTe 

Contact arrière 


Fig.II.ll : Cellule solaire au tellurure de cadmium. 


Malgré des rendements excédant les 21% sur une surface de 1 cm 2 . Cependant les problèmes 
d’environnement associés à l’utilisation du cadium freinent les tentatives de développement de 
cette filière [11]. 

II. 10.3 Filière Cu(In,Ga)Se 2 

L’émergence des cellules solaires, basées sur les matériaux à structures chalcopyrites, a fait 
suite à la réalisation de détecteurs photovoltaïques CuInSe 2 / CdS par S. Wagner et Coll de Bell 
Téléphone en 1974-1975. S. Wagner et Coll avaient rapporté la préparation d’hétérojonctions 
p-n CuInSe 2 /CdS à partir d’un monocristal CuInSe 2 , dont la conductivité est de type p, sur 
lequel une couche mince de CdS de 5-10 pm d’épaisseur avait été déposée. La réponse 
photo voltaïque de ces détecteurs pour une lumière incidente à travers la fenêtre en CdS donnait 
lieu à de très hauts rendements quantiques (70%) uniformes entre 0,55 et 1,25 pm. Ces 
hétéroj onctions avaient des rendements solaires de l'ordre de 5%. Les composés ternaires 
chalcopyrites qui peuvent jouer le rôle d’absorbeur sont principalement CuGaSe 2 , CuInS 2 , 
Cu(In, Al) Se 2 et Cu(In, Ga)Se 2 . Du fait de leurs grands coefficients d’absorption optique, une 
épaisseur de 1,5-2 pm est suffisante pour absorber la partie utile du spectre solaire [19]. 
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La Figure 11.12 représente les caractéristiques I(V) réelle et simulée d’une cellule CIGS 
réalisée en laboratoire avec un rendement de 17.5% pour un ensoleillement AMI. 5. 
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Fig.II .12 : Caractéristique courant-tension d’une cellule CIGS [18], 

Bien que la majeure partie de la production de cellules solaires soit faite sur des substrats de 
verre, on voit apparaître progressivement de nouveaux substrats souples compatibles avec la 
filière couche mince. Ces nouveaux substrats (métal ou plastique) présentent les avantages 
d’être plus légers que le verre, plus faciles à mettre en œuvre et surtout moins coûteux. Ces 
caractéristiques font d’eux des substrats de choix pour le futur [18]. 

II.10.4 Multi-jonctions en couches minces 

Le principe d’une cellule solaire multi-jonction étant d’absorber plus de longueurs d’onde 
du spectre solaire, naturellement l’idée viendrait d’ajouter des diodes supplémentaires avec des 
matériaux à gaps décroissants pour augmenter encore plus le rendement [20] (fig.II. 13). 
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Egi > E g 2 > E g 3 

Cellule 1 (E g i) 
Cellule 2 (E g2 ) 
Cellule 3 (E g j) 




Fig.II.13 : Principe de la cellule solaire multi-j onctions [18]. 


Pour les cellules multi-jonctions en semi-conducteurs III-V (GaAs, AlGaAs, InGaAsP), de 
très bons rendements sont obtenus avec cette filière mais les coûts de fabrication sont plus 
élevés (bâti sous ultravide, croissance très lente). Ces cellules sont destinées à des applications 
spatiales, un rendement de 32 % a été obtenu avec une triple jonction GalnP/GaAs/Ge sans 
concentrateur et 41 % avec concentrateur [9]. 

11.11 Cellule solaire en CIGS 

Le développement de cellules solaires à base de CIGS a débuté dans les années 1970. Les 
rendements de conversion des dispositifs dépassent 22% en laboratoire et atteignent 16% à 
l’échelle industrielle [21], Bien qu’une cellule CIGS soit constituée d’au moins 6 matériaux 
différents, c’est principalement la compréhension et l’amélioration des propriétés du matériau 
absorbeur qui sont au cœur des progrès réalisés. L’utilisation de substrats métalliques à la place 
du verre constitue une véritable problématique additionnelle. 

11.12 Structure de la cellule CIGS 

Dans sa structure la plus répandue, une cellule CIGS est formée d’un empilement de 
plusieurs matériaux en couches minces déposés successivement sur un substrat. La figure 
(11.14) présente la structure standard d’une cellule à base de Cu(In, Ga)Se 2 . 

Les principaux éléments dans la cellule sont [5] : 
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II.12.1 Substrat 

Le substrat le plus couramment utilisé est le verre sodo-calcique, qui présente toutes les 
propriétés requises. Il contient notamment du sodium, qui diffuse vers le CIGS lors de l’étape 
de recuit. Son effet bénéfique sur le fonctionnement de la cellule a été démontré par de 
nombreux auteurs [22,23]. 



(a) (b) 


Fig.II.14 : (a) schéma de l’empilement des différentes couches d’une cellule en CIGS, (b) 
image de la tranche d’une cellule en CIGS obtenue par un microscope électronique à balayage 
(MEB). 

II. 12.2 Contact arrière 

Le contact arrière, qui constitue l’anode de la cellule solaire, est un métal déposé sur la partie 
inférieure de l’absorbeur dont le rôle est de collecter les porteurs à partir de l’absorbeur et de 
les livrer à la charge externe. Il doit donc, répondre à certains critères chimiques, mécaniques 
et économiques. De nombreux métaux ont été étudiés comme le Tungstène, le Molybdène, le 
Chrome, le Tantale, le Manganèse [24], l’Or, l’Aluminium, l’Argent et le cuivre [25], Parmi 
ces métaux, le Molybdène qui est le plus couramment utilisé comme contact arrière des cellules 
solaires CIGS par les laboratoires de recherche et en industrie, en raison de sa stabilité relative 
aux températures élevées rencontrées lors du procédé de fabrication du CIGS (de Tordre de 
550 à 600°C), et de sa faible résistance de contact avec le CIGS. Son épaisseur est d’environ 
300 nm à 1000 nm. La méthode de dépôt utilisée est la pulvérisation cathodique. La valeur 
typique de la résistivité de Mo est de Tordre 5><10" 5 G. cm ou moins. De plus, il ne réagit pas 
avec le CIGS ou y diffuse. 
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En effet, le Mo peut réagir avec le Se lors du dépôt de CIGS pour former une couche de 
quelques nanomètres du MoSe 2 . Par conséquent, la structure CIGS/Mo devient alors 
CIGS/MoSe 2 /Mo. Le MoSe 2 est un semi-conducteur possédant un gap de 1,41 eV [26], il est 
plus élevé que celui du CIGS est permet ainsi, de diminuer les recombinaisons à l’interface de 
ces deux couches [27]. De plus, il a pour effet de donner un comportement ohmique à l’hétéro 
contact CIGS/Mo [26, 28], 

11.12.3 Absorbeur en CIGS 

L’absorbeur est la partie de la cellule dans laquelle les photons vont être convertis en paires 
électrons-trous. Dans le cas présent, il s’agit d’un semi-conducteur de type P de la famille I- 
III- VI, dont les différentes variantes déjà étudiées conduisent à la formule générale suivante : 
(Cu, Ag) (In, Ga, Al) (Se, S, Te) 2 . La substitution du Cu par de l’argent, de l’indium par du 
gallium ou de l’aluminium, du sélénium par du soufre ou du tellure, permet de modifier les 
propriétés du matériau comme sa largeur de bande interdite. Le composé le plus développé en 
recherche et en industrie, et qui constitue le sujet de notre étude, est le Cu(In, Ga)Se 2 [29]. 

11.12.4 Couche tampon en CdS 

La couche tampon est un semi-conducteur ayant une conductivité de type n dont le gap doit 
être supérieur à celui de l’absorbeur. Deux rôles lui sont majoritairement attribués : un rôle 
électrique, et un rôle de couche protectrice. Du point de vue électrique, cette couche adapte la 
largeur de la bande interdite entre l’ absorbeur et la couche fenêtre [30] et limite les 
recombinaisons des porteurs à l’interface de ces deux couches. Elle permet également de 
protéger la surface de l’ absorbeur lors du dépôt par pulvérisation cathodique de la couche de 
ZnO, qui peut engendrer des défauts à la surface du CIGS [31], 

Plusieurs semi-conducteurs de type n dont les bandes interdites sont entre 2,0 et 3,6 eV ont 
été appliquées comme tampon pour les cellules solaires CIGS. Cependant, le CdS déposé par 
bain chimique (CBD) reste préféré car il donne toujours des hauts rendements [32- 34], 
Néanmoins, sa largeur de bande interdite, égale à 2 ,4 eV relativement faible engendre des 
pertes optiques, de sorte que des quantités importantes de lumière sont absorbées avant de 
pouvoir atteindre la couche CIGS, même pour les couches de CdS aussi minces que 50 nm. 

Les photons absorbés dans le CdS ne contribuent pas au photocourant, à cause de la petite 
longueur de diffusion des porteurs minoritaires (trous) dans le CdS, ce qui conduit à un taux 
élevé de recombinaison. Pour cette raison, et en plus des contraintes liées à l’utilisation du 
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Cadmium (Cd) qui est un élément toxique, de nombreuses études sont menées pour substituer 
le CdS par d’autres matériaux, parmi eux, on peut citer (OH, S), Zn (O, OH), ZnS [31, 35] ou 
ZnSe [35-37], 

II. 12.5 Couche fenêtre en ZnO 

La couche fenêtre doit être à la fois transparente, afin de laisser suffisamment de lumière à 
la région de jonction et la couche absorbante du dispositif, et conductrice pour permettre le 
transport du photocourant généré au circuit extérieur sans trop de pertes. Cette couche fenêtre 
est constituée d’oxyde transparent conducteur (TCO), sur lequel une fine grille métallique 
doublée (Ni- Al) est déposée de façon à réduire la résistance série de la couche fenêtre et de 
collecter les charges générées par le dispositif. Le Ni sert de couche d’accroche et permet 
d’éviter l’oxydation de l’Al liée à la présence sous-jacente du TCO. Le TCO, généralement 
utilisé pour les cellules CIGS, est l’oxyde de Zinc (ZnO) déposé par pulvérisation cathodique, 
mais l’oxyde d’indium dopé à l’étain (ITO) peut être également rencontré [38], 

Une bicouche ZnO est souvent utilisée pour le contact frontal transparent. La première 
couche est une couche mince de ZnO intrinsèque (i-ZnO), qui est typiquement de (50 -100 nm) 
d’épaisseur et très résistive. Elle sert à limiter la formation de court-circuit dans les zones 
présentant un recouvrement imparfait du CIGS par la couche tampon [39]. La seconde couche 
est fréquemment du ZnO dopé à l’Aluminium (ZnO : Al), mais d’autres dopants comme le 
Bore (ZnO : B) [40] ou le Gallium (ZnO : Ga) [41] existent aussi. L’épaisseur et le dopage de 
cette couche sont ajustés de façon à combiner transparence et conductivité. La bande interdite 
de la couche fenêtre doit être aussi élevée et aussi mince que possible afin de maintenir une 
faible résistance série. Le ZnO et le ZnO : Al ont respectivement des bandes d’environ 3,2 eV 
et 3,6 à 3,8 eV. 
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II.12.6 Jonction p-n et diagramme d’énergie 

La jonction p-n dans les cellules solaires CIGS est formée entre le semi-conducteur type p 
(CIGS) et le type n (CdS). Le diagramme d’énergie de bande de cette cellule est représenté sur 
la figure (11.15). La différence des affinités électroniques des deux matériaux constituant la 
jonction se traduit par une discontinuité de bande interdite. La courbure de la bande de 
conduction ou ÀEc (Conduction Band Offset) entre le CIGS et le CdS est par convention 
positive, ce qui signifie que la bande de conduction de CdS est supérieure à celle de CIGS. 



02 04 2 4 2 5 2 6 2 7 

Distance [pm] 


Fig.II.15 : Structure de bande d’une cellule solaire Cu(In, Ga)Se 2 . 

La valeur optimale de la courbure de la bande de conduction est AEc = 0,3 eV est obtenue 
pour un gap de CIGS=1,15 eV. Une courbure de la bande de conduction supérieure à 0,5 eV 
crée une barrière forte pour les électrons photo- générés limitant ainsi, le transport du courant. 
En revanche, une courbure de la bande de conduction nulle ou négative conduit à une 
augmentation de la recombinaison à l’interface, ce qui limite la tension de la cellule. 

La présence d’une couche de type inversé (type n) identifiée comme étant une couche de 
défauts de surface ou SDL (Surface Defect Layer) à la surface du matériau CIGS était rapportée 
par plusieurs auteurs [42], Elle possède un gap plus large que le CIGS, ce qui permet de réduire 
la recombinaison à l’interface CdS/CIGS et contribue ainsi à l’amélioration des performances 
de la cellule [43, 44]. 
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Le tableau II. 1 regroupe certains paramètres physiques de la cellule à base de CIGS. 


Propriété 

CIGS 

CdS 

ZnO 

E g (eV) 

Varie (1.011/1.68) 

2.4 

3.3 

Xe (cm.- 3 ) 

4.58 

4.4 

4.4 

£ r 

13.6 

10 

9 

N c (cm~ 3 ) 

2.2xl0 18 

2.2xl0 18 

2.2xl0 18 

N v (cm~ 3 ) 

1.8xl0 19 

1.8xl0 19 

1.8xl0 19 

p n (cm 2 /Vs) 

100 

100 

100 

g p (cm 2 /Vs) 

25 

25 

25 

N a Ou N d (cmT 3 ) 

N a =2xl0 16 

N d =10 17 

10 17 

m e 

0.2m 0 

0.9m 0 

0.9m 0 

m t 

0.8m 0 

0.8m 0 

0.8m 0 

Tn OU Tp (s) 

T n =2xl0 9 

T p =10 13 

/ 

S n ou S p (cm/s) 

S n =10 7 

S p =10 7 

/ 


Tableau II. 1 : Paramètres physiques des matériaux de la cellule solaire en CIGS [45]. 

11.13 Contenu de Ga dans la cellule solaire en CIGS 

Le rendement d'une cellule solaire est déterminé à partir de la caractéristique I-V de la 
cellule solaire ainsi la tension en circuit ouvert, le courant de court-circuit et le facteur de forme 
qui constituent les facteurs de qualité de la cellule solaire. Une augmentation de la bande 
interdite de l’absorbeur en augmentant ainsi le contenu de Ga dans x=Ga/(In+Ga), conduit à 
une augmentation proportionnelle de V co . L’accroissement de la bande interdite entraîne une 
augmentation de la tension en circuit ouvert et du facteur de forme. 

La réponse spectrale d'une cellule solaire en CIGS avec la variation de la composition Ga 
est montrée sur la figure 11.16. La réduction des performances avec une augmentation de la 
composition de Ga peut se justifier du fait de l’augmentation de la bande interdite de 
l’absorbeur et il y’aura moins de la contribution de photons, qui pénètrent plus profondément 
et générant le photocourant. 
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Fig.II.16 : Réponse spectrale d’une cellule en CIGS pour différentes proportions x de Ga 
[46], 


Des exemples de caractéristiques I-V pour trois bandes interdites différentes sont 
représentés sur la figure 11.17. On remarque que plus la bande interdite augmente, la tension en 
circuit ouvert augmente aussi contrairement au courant de court-circuit qui diminue. 



Fig. 11.17 : caractéristiques I-V pour trois bandes interdites différentes [47]. 
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11.14 Avantages de la fabrication des cellules solaires en CIGS 

L'utilisation de cellules solaires en CIGS a le potentiel de diriger le marché des cellules 
solaires en couches minces. Cela est dû aux nombreux avantages que les cellules en CIGS 
offrent à la fois au fabricant et au consommateur. Comme la recherche continue à améliorer à 
la fois le rendement des cellules et la réduction des coûts de fabrication, les cellules solaires en 
CIGS ont le potentiel d'acquérir la majorité de la part du marché des couches minces [48]. 

II. 14.1 Fabrication peu coûteuse 

Un des principaux avantages de cellules en CIGS est qu’elles peuvent être commercialement 
fabriquées à une fraction du coût par rapport aux cellules solaires en Si. Les cellules cristallines 
de Si doivent être fabriquées avec le silicium extrêmement pur. Une fois qu’un fabricant a 
abouti au silicium avec une assez grande pureté, le matériau doit être chauffé à des températures 
extrêmement élevées pour se former dans une structure cristalline. Le lingot du silicium 
cristallin est ensuite découpé en plaquettes qui doivent subir de nombreuses étapes de gravure 
et de dopage avant qu’un produit final soit atteint [49]. Ce processus est très long et coûteux 
pour le fabricant. 

Les cellules solaires en couches minces, comme en CIGS peuvent être produites de manière 
beaucoup plus efficace. Le matériau CIGS peut être déposé sur le substrat et le contact arrière 
par un procédé appelé déposition qui permet aux éléments de s’adhérer sur le substrat à une 
température suffisamment élevée. Ce processus ne nécessite pas des tolérances élevées par 
rapport aux cellules en Si et il est globalement un processus plus facile pour la fabrication [50]. 
En plus des avantages de la déposition, le CIGS peut être fabriqué avec une technique appelée 
sur rouleaux (" roll-to-roll "). La fabrication " roll-to-roll " est illustrée sur la figure 11.18. 

Avec la technique de fabrication "roll-to-roll ", les sociétés peuvent produire des cellules 
solaires à grande échelle. Une fois le substrat flexible est déroulé, il subit plusieurs processus 
qui consistent à déposer l’ensemble des couches sur le substrat. Le processus réduit 
considérablement le temps de production mais permet également un cycle de production 
continu. Plutôt que de couper des plaquettes individuelles de Si, la fabrication roll-to-roll 
supporte une feuille continue de cellules solaires produites. Ce processus est si efficace qu'il a 
permis à des sociétés comme Miasolé et Solibro d'offrir des cellules solaires CIGS avec un prix 
de 1 $ / W [51]. 
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Fig.II.18 : Le processus de fabrication " roll-to-roll " de la cellule solaire CIGS. 

II. 14.2 Flexible 

La flexibilité des cellules CIGS produites sur un substrat métallique mince est aussi un 
avantage énorme. Cette flexibilité donne au CIGS un avantage par rapport aux cellules Si 
conventionnelles en ce qui concerne le montage. Une cellule Si standard possède des propriétés 
compatibles avec un panneau en verre qui rend une cellule Si rigide mais cassante, ceci limite 
les options de mise en œuvre de panneaux solaires. La nature flexible de cellules solaires CIGS 
leur permet d'être montées sur des structures qui ne peuvent pas être parfaitement plates. La 
recherche a montré que le montage des cellules CIGS est très pratique et avec une grande dureté 
sur des avions sans pilotes [52], 

II. 14.3 Poids léger 

Le sous-produit d'être mince et flexible est aussi léger. Étant léger a ses avantages pour les 
personnes qui souhaitent prendre leurs cellules solaires avec eux lorsqu’ils voyagent. Leur 
portabilité reste hors de portée pour les cellules à base de Si. 

II. 14.4 Etat actuel de la technologie 

La recherche scientifique continue d'améliorer le rendement observé dans les cellules en 
CIGS fabriquées. Le rendement record du laboratoire national d'énergies renouvelables 
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(NREL) pour de nombreuses technologies de cellules solaires par rapport à la date à laquelle 
le rendement a été enregistré est présenté dans la figure 11.19. 
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Fig.II.19 : Meilleurs rendements des cellules solaires du laboratoire NREL en fonction du 
temps [53], 


11.15 Conclusion 

Dans le présent chapitre, nous avons présenté les bases fondamentales à la compréhension 
de la cellule solaire. Nous avons rappelé quelques notions sur le rayonnement solaire, et son 
application dans le domaine photovoltaïque. Nous avons ensuite expliqué le fonctionnement 
des cellules photovoltaïques et leurs caractéristiques principales ainsi que les paramètres qui 
influent sur leurs caractéristiques. Nous avons aussi présenté la structure et quelques notions 
sur les cellules solaires à base de CIGS. 
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III. 1 Introduction 

SILVACO (Silicon Valley Corporation) est une société Américaine, ayant son siège à Santa 
Clara en Californie. Elle est un des principaux fournisseurs de chaînes professionnelles de 
logiciels de simulation par éléments finis et de conception assistée par ordinateur pour les 
technologies de l’électronique TC AD (Technology Computer Aided Design). Ces outils sont 
employés par les compagnies de microélectronique dans le domaine de la recherche, du 
développement et de la conception de dispositifs [1]. 

Les simulateurs TCAD permettent la modélisation du comportement électrique d’un 
composant électronique, dans le but d’économiser le temps et le coût de développement, et 
ainsi de pouvoir envisager et optimiser des solutions pour améliorer les performances des 
dispositifs. Les simulateurs existants sur le marché (ISE, SYNOPSIS, SILVACO, ...) ont le 
même principe de fonctionnement [2]. 

Dans ce chapitre, nous allons présenter le logiciel de simulation TCAD-SILVACO, son 
principe de fonctionnement, et ces modules de simulation : ATHENA pour la simulation des 
processus technologiques [3], et ATLAS pour la simulation du comportement électrique des 
composants [4], Nous allons aussi donner les différentes commandes (instructions) et les bases 
physiques du principe de fonctionnement du simulateur ATLAS. 

III.2 Présentation du package des programmes SILVACO 

SILVACO présente un ensemble des outils de simulation interactifs permettant la 
conception et l’analyse de la plupart des dispositifs à semiconducteurs VWF (Virtual Wafer 
Fab) [1]. Les composants de base de VWF (Figure III. 1) sont [1,5] : 

1. Les outils de simulation (VWF CORE TOOLS). Ces outils simulent soit leurs processus 
de fabrication ou soit leurs comportements électriques. Les outils de simulation sont ATHENA, 
ATLAS et SDUPEM3. 

2. Les outils interactifs (VWF INTERAVTIVE TOOLS). Ces outils sont désignés pour 
être utilisés en mode interactif dans la construction d’un seul fichier d’entrée. En étant basé sur 
une interface utilisateur qui est graphique (Graphical User Interface, GUI). Ainsi le travail de 
construction du fichier d’entrée devient plus efficace. Les outils interactifs peuvent être utilisés 
soit en relation avec un ensemble de fichiers, ou comme des composants intégrés dans 
l’environnement des outils d’automatisation. 
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3. Les outils d’automatisation (VWF AUTOMATION TOOLS). Ces outils permettent à 
l’utilisateur d’exécuter sur une grande échelle des études expérimentales afin de créer des 
résultats pour l’analyse statistique suivante. Les outils automatiques se servent de la 
technologie de base de données réparties et des méthodes de logiciels de développement 
d’interprocessus. 



Fig.III.l : Organigramme de la structure VWF [1], 
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111.2.1 Les outils de simulation 

> ATLAS : simulateur physique de dispositifs à semiconducteur 2D ou 3D qui permet de 
simuler le comportement électrique de structures spécifiées des dispositifs à 
semiconducteur [1], 

> ATHENA : simulateur 2D de procédés technologiques qui permet de développer et 
optimiser des processus de fabrication de semiconducteur (les différentes étapes effectuées 
en Salle Blanche). Il fournit une plate-forme pour simuler l'implantation ionique, la 
diffusion, la gravure, le dépôt, la lithographie et l'oxydation. Il remplace des expériences 
coûteuses par des simulations [1]. 

> SSUPREM3 : simulateur de procédé 1D permet la prédiction des profils de dopage et les 
épaisseurs produits par le traitement des semi-conducteurs, il est capable de simuler un 
flux complet d'étapes de processus [1], 

111.2.2 Les outils interactifs 

> DECKBUILD : environnement où est défini le programme de simulation. DeckBuild 
offre une aide et une prise en charge qui permet d’éviter l’apprentissage d’une syntaxe de 
commande généralement très complexe. Il permet aussi une exécution automatique et très 
pratique des programmes et des sauvegardes systématiques des étapes intermédiaires 
d’exécution [1], 

> TONYPLOT : environnement où sont visualisés les résultats des simulations (structure 
du composant, distribution de grandeurs divers dans celui-ci, caractéristiques 
électriques...). 

> DEVEDIT : environnement où est dessinée la structure. On peut créer des nouvelles 
structures ou même modifier des structures existantes. On peut aussi définir des maillages 
ou raffiner les maillages existants [1], 

> MASKVIEWS : outil de dessin des masques (éditeur de layout) [1]. 

> OPTIMIZER : optimise les paramètres de la structure de façon à obtenir en final la valeur 
du paramètre que nous lui avons définie au préalable [1], 

> MANAGER : outil de gestion des fichiers utilisés et créés par VWF [1], 

> SPDB : (Semiconductor Process Data Base), c’est un produit séparé, il n’est pas un outil 
interactif, mais il peut être utilisé avec DeckBuild. Il a été conçu pour stocker un grand 
nombre de profils de dopage mesurés expérimentalement ainsi que les données qui 
décrivent les conditions des expériences [1]. 
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III.3 Présentation d’ ATLAS 

Le logiciel de simulation ATLAS (de SILVACO) permet la modélisation 2D et 3D des 
composants. Cet outil sert à la modélisation des dispositifs à semi -conducteur avant leur 
fabrication. Il est très utile dans le développement de beaucoup de projets de recherches. 

Dans la figure III.2, nous présentons les informations d’entrée et de sortie pour ATLAS. 


DevEdit 



Athéna 



Fichier de structure 


DeckBuild — ► Fichier de commande 


I 



Entrées 


Sorties 


Fig.III.2 : Entrées/sorties du simulateur Atlas [6], 


La plupart des simulations réalisées par ATLAS utilisent deux fichiers d’entrée [6] : 

> Un fichier script contenant les commandes pour que ATLAS s’exécute (environnement 
DECKBUILD utilisé), 

> Un fichier structure définissant l’empilement à simuler. 

Le simulateur fournit 3 types de fichiers de sortie : 

> Un fichier Runtime qui donne la progression, les erreurs et les messages d’avertissements 
pendant la simulation. 

> Un fichier log stockant toutes les valeurs de tensions et des courants provenant de l’analyse 
du dispositif simule (c’est le fichier du comportement électrique), 

> Un fichier solution stockant les données 2D et 3D correspondant aux valeurs de variables 
spécifiques à l’intérieur du dispositif, à une tension donnée. 
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Les deux derniers fichiers sont traités par l’outil de visualisation TonyPlot. Tonyplot permet 
de visualiser la structure en 2D ou 3D et les caractéristiques électriques demandées. 

III.4 Commandes dans un programme ATLAS 

Nous allons maintenant présenter l’ordre des commandes fondamentales dans un 
programme ATLAS (Figure III.3). Ainsi il existe cinq groupes de commandes, ces groupes 
doivent être organisés correctement. Si l’ordre n’est pas respecté, un message d’erreur apparaît 
et le programme ne s’exécute pas d’une façon correcte. 


Groupe 


Commandes 


1 . Spécification de la structure 


> 


2. Spécification des modèles de couches ^ 


3. Sélection de la méthode numérique 


4. Spécification des solutions ^ 


5. Analyses des résultats 
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Fig.III.3 : Les commandes fondamentales dans un programme ATLAS [1], 
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III.4.1 Spécification de la structure 

Il s’agit de définir la structure la plus réaliste. La construction de la structure consiste en 
quatre étapes : définition du maillage, définition des régions, définition des électrodes et 
définition des dopages, d'après [7,8] : 

> MESH (définition du maillage) : Le maillage divise la structure simulée en petites cellules 
pour pouvoir résoudre les équations fondamentales numériquement. L’élément de maille 
utilisé est le triangle. Pour assurer une bonne précision lors de la simulation, il faut avoir 
le maillage le plus fin possible. Par contre si le maillage est très fin, le temps de calcul est 
très élevé, car il y a plus d’éléments à calculer. Le format général pour définir la maille 
est : 


X. MESH LOCATION = <VALEUR> SPACING = <VALEUR> 

Y. MESH LOCATION = <VALEUR> SPACING = <VALEUR> 

La figure III.4 montre le maillage de notre structure. 
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Fig.III.4 : Maillage du dispositif. 
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> REGION (définition des régions) : Après avoir défini le maillage, il est nécessaire de 
définir les régions dont le format de définition des régions avec dimensions en micromètre 
est le suivant : 

REGION nombre = <integer><material_type>< position des paramètres > 

La figure III.5 montre la façon de définition des matériaux pour chaque région. Notons 
que le code de couleur indique le matériau. Les régions ont des lignes horizontales pour 
marquer leurs limites. 
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Fig.III.5 : Définition des régions et des matériaux par le logiciel ATLAS. 


> ELECTRODE : Atlas a une limite de 50 électrodes qui peuvent être définies. Le format 
de définition d’électrodes est co mm e suite : 

ELECTRODE name = < nom de l’électrode >< position des paramètres > 
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Les positions des électrodes sont localisées par les instructions suivantes : 

X.MIN : Spécifie le point de départ de l’électrode. 

RIGHT : la position de l’électrode se situe à la droite de la structure (inverse : LEFT). 
TOP : la position de l’électrode est en haut de la structure (inverse : BOTTOM). 

La figure III.6 montre la position de l’anode et la cathode dans la structure. 
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Fig.III.6 : Définition des électrodes par le logiciel ATLAS. 

> DOPING (définition des dopages) : Le dernier aspect de la construction de la structure est 
le dopage. Le dopage peut être de type n ou p, aussi la distribution peut être uniforme, 
gaussien, etc. . . 

Le format de la déclaration de dopage dans ATLAS est comme suite : 

DOPING< type de distribution >< type de dopant >< position des paramètres > 
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Fig.III.7 : Distribution du dopage dans une structure par le logiciel ATLAS. 

III.4.2 Spécification des modèles 

Après la construction de la structure, il est nécessaire de préciser les modèles physiques et 
de définir les matériaux. La spécification de modèles et matériaux comporte les étapes [5, 9, 
10 ]: 

> MATERIAL : Le format de la déclaration du matériau est le suivant : 

MATERIAL < localisation >< définition du matériau > 


Exemple : 

MATERIAL MATERIAL=CdS EG300=2.42 PERMITTIYITY =10 \ 
AFFINITY=4.5 MUN=100 MUP=25 NC300=2.2el8 NV300=1.8el9. 

Plusieurs autres paramètres peuvent être définis par la déclaration MATERIAL, comme 
exemples de ces paramètres : l’énergie du gap (EG300), la permittivité (PERMITIVITY), 
l’affinité (AFFINITY), la mobilité d'électrons (MUN), la mobilité de trous (MUP), la 
densité d’états en Ec (NC300) et la densité d’états en Ev (NV300). 
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> MODELS : Ils indiquent les inclusions de différents mécanismes physiques 
(recombinaison, mobilité, etc...)- Le choix du modèle dépend des matériaux choisis pour 
la simulation. 

La syntaxe de la déclaration du modèle est la suivante : 

MODEL< paramètres générales >< paramètres du modèle > 

Exemple : S RH CONMOB FLDMOB. 

SRH est le modèle de Shockley-Read-Hall. CONMOB est le modèle de la mobilité 
dépendante de la concentration du dopage. FLDMOB est le modèle de la mobilité 
dépendante du champ électrique. 

> CONTACT : Il existe aussi les modèles qui concernent le type de contact soit contact 
ohmique qui connaît par défaut et le contact redresseur (contact de Schottky), ce contact 
est reconnu par un travail de sortie (work function) comme : 

Contact name= < integer > work=< integer > 

Voici un exemple des contacts déclaration : 

Contact name =anode workfun= 4.74 

> INTERFACE : Les limites des semiconducteurs ou des isolants sont déterminées par la 
déclaration « interface ». La syntaxe est la suivante : 

INTERFACE <paramètres> <interface_dimensions > 

L’exemple suivant montre l’utilisation de la déclaration d’interface : 

Interface x.min=0 x.max=2000 y.min=0 y.max=0 s.n=le3 s.p=le3 

Les valeurs max et min déterminent les limites, s.n et s.p sont les vitesses de 
recombinaison superficielle des électrons et des trous aux interfaces. 


> LUMIERE : la lumière est modélisée dans ATLAS par la déclaration : 

BEAM <parameters> 
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En voici un exemple : 

beam num=l x.origin =1000 y.origin=-30 angle=90.0 power. file=AMOdis.spec reflects=l 

front.refl back.refl 

Cette déclaration définit le nombre des faisceaux lumineux qui est 1, l’origine du rayon 
pour des sources optiques 2D ; le x.origin et le y.origin sont les coordonnées x et y de 
l’origine du faisceau lumineux en microns. L’angle d’incidence de 90 degrés indique que 
le rayon est perpendiculaire à la surface éclairée de la cellule solaire, power .file spécifie 
la source spectrale de la lumière en indiquant le nom d’un fichier d’extension . spec qui 
contient les valeurs des intensités du spectre utilisé en fonction de leurs longueurs d’onde. 

Le paramètre reflects=l donne le nombre de réflexions du faisceau lumineux. La 
déclaration front.refl et back.refl permet au rayon lumineux de se réfléchir en face avant 
et arrière de la cellule. 

III.4.3 Sélection de la méthode numérique 

Après la précision du modèle de matériaux, la sélection de la méthode numérique doit être 
indiquée. Parmi les différentes méthodes numériques de résolution des systèmes d’équation, 
trois types de techniques sont utilisée dans « Silvaco- Atlas » : 

> Méthode de Gummel 

> Méthode de Newton 

> Méthode des blocks 

Pour la méthode de GUMMEL, elle est utilisée pour la résolution des équations a une seule 
inconnue et des variables fixées constantes, la résolution se faite d’une manière itérative 
jusqu’à la solution est achevée. La deuxième méthode de NEWTON est la méthode la plus 
répandu dans le programme Atlas, elle est utilisée pour résoudre des systèmes d’équations à 
plusieurs inconnues. La troisième méthode est une combinaison des deux méthodes, elle est 
utiles pour les systèmes d’équation mixtes. 
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III.4.4 Spécification des solutions 

Après avoir sélectionné la méthode numérique il est nécessaire de déterminer les solutions. 
La spécification de solution est donnée par les déclarations : LOG, SOLVE, LO AD, et SAVE, 
d’après [7,9] : 

> Log : L’instruction Log est utilisée pour enregistrer toutes les caractéristiques terminales 
dans un fichier de données. L'exemple suivant montre un exemple d’utilisation de 
l’instruction LOG. 

logoutfile = myoutputfile.log 

L'exemple montre l’enregistrement de données dans le fichier (myoutputfile.log). 

> Solve (Résoudre) : L’instruction SOLVE suit l'instruction LOG. Elle effectue une solution 
pour un ou plusieurs points de polarisation. Cette déclaration a le format de base suivant : 

SOLVE <bias_points> 

Par exemple : 

solve vanode=0 vstep=0.01 vfinal=2 name=anode 

> Load et Save : La commande LOAD (chargement) introduit, à partir d’un fichier 
existant, pour chaque point de polarisation les solutions précédentes en tant que 
proposition initiale. L’instruction SAVE (sauvegarde) permet l’enregistrement de toutes 
les informations obtenues pour un nœud dans un fichier de sortie. Les instructions pour 
load et save sont comme suit : 


LOAD INFILE=<filename> 
SAVE OUTFILE =<filename> 


III.4.5 Analyse des résultats 

Une fois la solution a été trouvée pour un problème de dispositif à semiconducteurs, les 
informations peuvent être affichées graphiquement. D’après [1] : 

> EXTRACT : les commandes de ce type sont utilisées pour extraire les valeurs bien 
précises des paramètres des deux types des fichiers log ou structure. 

> TONYPLOT : démarre le programme « TonyPlot » permettant une représentation 
graphique des données. 
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III.5 Les équations de Base de Semi-conducteur en ATLAS-SILVACO 

Ce simulateur de dispositifs est basé sur la résolution simultanée de l’équation de Poisson 
et des équations de continuité pour les électrons et les trous. Il calcule à chaque instant et en 
tout point de l’espace en une suite d’éléments finis, la concentration des électrons et des trous 
et la valeur du potentiel électrostatique [8]. 

L’équation de Poisson montre une relation entre le potentiel et la densité des porteurs [4] : 

A V = -^[p -n + N£ - A/f] III. 1 

Où V le potentiel, q la charge élémentaire d’électrons, £ la constante diélectrique, et 
les concentrations des dopants donneurs et accepteurs ionisées, et n et p les densités des 
électrons et des trous. 

Les équations de continuité pour les électrons et les trous sont [4] : 
dn 1 . ^ 

— — = H — divj n + G n — R n III. 2 

dt q 

dp 1 

— = — divj v + G v — R v III. 3 

ot q F w y 

Où J n et J p sont les densités de courant des électrons et des trous, G n et G p les taux de 
génération pour les électrons et les trous, R n et R p les taux de recombinaisons pour les électrons 
et les trous. 

Les dérivations basées sur la théorie de transport de Boltzmann ont prouvé que la densité 
des courants dans les équations de continuité peut être définie par le modèle d’entrainement 
diffusion. 

Les densités de courants sont exprimées en deux termes. Le premier correspond à 
l’entrainement des charges par le champ électrique, le second correspond à la diffusion des 
porteurs [4]. 

J n = qnp n E + qD n grad(n ) III. 4 

III. 5 
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Où D n et D p sont les coefficients de diffusion. Ils sont directement reliés aux mobilités des 
porteurs par la relation d’Einstein : 


kT 

D n = — I Un III. 7 

q 


_ kT 

Dp 


III. 8 


III.6 Conclusion 

Ce chapitre a présenté quelques informations de base en programmation par Silvaco-Atlas, 
alors plusieurs aspects spécifiques à la méthodologie de programmation ont été présentés 
comprenant structure de programme, Spécification des paramètres des matériaux et modèles 
physiques et sélection de la méthode numérique. 
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IV. 1 Introduction 

Dans ce chapitre nous présentons les résultats de l’étude par simulation numérique des 
caractéristiques électriques d’une cellule solaire en couches minces et à hétéroj onction à base 
de CIGS et leurs interprétations. En premier lieu, nous allons simuler deux cellules solaires de 
structures : ZnO/CdS(n)/CIGS(p)/Mo et ZnO/ZnS(n)/CIGS(p)/Mo en utilisant le logiciel 
Silvaco Atlas-2D. Ensuite nous allons comparer les résultats de la simulation électrique de la 
cellule en CdS avec ceux trouvés expérimentalement. 

En deuxième lieu, nous allons comparer les résultats de la simulation électrique de la cellule 
en CdS avec la cellule en ZnS. Après cela, nous choisirons la meilleure cellule pour étudier 
l’effet du dopage et de l’épaisseur des couches tampon /absorbeur sur les paramètres 
photovoltaïques (J cc , V co , FF, t] ) de la cellule solaire dans le but d’aboutir à un rendement 
maximal de la cellule solaire à base du CIGS. 

IV.2 Cellule solaire en CdS/CIGS 

IV.2.1 Structure de la cellule en CdS/CIGS 

Nous avons pris comme référence la structure développée par Philip Jackson et al. [1], la 
cellule solaire en CIGS à étudier possède une structure p-n à hétérojonctions, la couche tampon 
en CdS de type (n) et la couche absorbeur en CIGS de type (p) avec deux électrodes l’anode 
en molybdène (Mo) et la cathode en ZnO. La construction de cette structure dans le simulateur 
Atlas consiste en quatre étapes : définitions du maillage, des régions, des électrodes et du 
dopage. 

=J- Première étape : nous avons défini le maillage de la structure, le maillage divise et 
discrétise la structure à simuler en petites cellules pour permettre la résolution 
numérique des équations fondamentales des semi-conducteurs (équations de 
continuité et de Poisson pour les électrons et pour les trous) afin de trouver la 
distribution des porteurs de charges, du potentiel électrique et du courant. Le maillage 
doit être le plus fin possible pour l’obtention de bons résultats de simulation. 

4- Deuxième étape : nous avons défini les régions avec ses matériaux représentant les 
couches formant la cellule solaire en ZnO/CdS(n)/CIGS(p)/Mo. 

i- Troisième étape : nous avons défini les régions et les matériaux spécifiés de chaque 
électrode. La cathode en ZnO est située sur la couche en CdS et l’anode en Mo est 
réalisée sur toute la surface inférieure de la couche en CIGS de la cellule solaire. 
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■4- Quatrième étape : nous avons défini le type de dopage et sa concentration, la couche 
en CdS de type n et la couche en CIGS de type p sont dopées uniformément avec des 
concentrations de 3xl0 17 cm" 3 et 8xl0 16 cm" 3 respectivement. 


La figure IV. 1 représente la structure simulée de la cellule solaire en CdS/CIGS. 


ATLAS 



0123456789 10 

Microns 


Fig.IV.l : Structure bidimensionnelle d’une cellule solaire en CdS/CIGS simulée 
dans Atlas-Silvaco. 
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IV.2.2 Propriétés des différents matériaux de la cellule solaire 

Nous avons utilisé les paramètres optiques et électroniques pour définir les différents 
matériaux qui composent la cellule solaire. 

Les paramètres des matériaux en ZnO, CdS et CIGS utilisés dans la simulation de la cellule 
solaire sont reportés dans la littérature [1-6] et sont groupés dans le tableau IV. 1. 


Régions 


Paramètres 


Epaisseur (pm) 


E g (eV) 



CIGS 



N c (cm 3 ) 


N v (cm 3 ) 

1.8xl0 19 

1.8xl0 19 

1.8xl0 19 

Mobilité des électrons 
(cm 2 /V.s) 

100 

100 

100 

Mobilité des trous 
(cm 2 /V.s) 

25 

25 

25 

Permittivité relative 

9 

10 

13.6 

Affinité 

4.7 

4.5 

4.8 


électronique (eV) 


Tableau IV.l : Paramètres de simulation la cellule solaire en CIGS. 


Les paramètres optiques ; la partie réelle et la partie imaginaire de l’indice de réfraction, 
n(À,) et k(X), du CdS et du CIGS se trouvent dans le logiciel Atlas tandis que celles du ZnO sont 
dans la référence [7]. 
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IV.2.3 Modélisation physique de la cellule solaire dans Atlas 

Les modèles que nous avons utilisés dans notre simulation sont décrits dans la 
documentation d’Atlas et dans la littérature [8-10]. Nous avons considéré le modèle de 
recombinaison Shockley-Read-Hall (SRH) pour les défauts. Le modèle de densité de défauts 
comprend deux distributions de défauts de formes gaussiennes pour les états donneurs pour 
ZnO [2], les états accepteurs correspondants au CdS [2] et les états donneurs pour CIGS [11]. 
Nous avons groupé les valeurs des paramètres de la densité de défauts où N DG ( AG ) est la densité 
des états donneurs (accepteurs) de formes gaussiennes, Ea(d) est l’énergie du pic de la 
gaussienne, W G la largeur de la distribution et a la section de capture (tableau IV.2). Le modèle 
de recombinaison à l’interface CdS/CIGS est aussi pris en considération dans notre simulation ; 
les vitesses de recombinaison des électrons (trous) 5 n (5 p ) = 10 5 cm/s. 



ZnO 

CdS 

CIGS 

N D g, N ag (cm 3 ) 

D :10 17 

A :10 15 

D :10 15 

E A .E D (eV) 

Milieu du gap 

Milieu du gap 

Milieu du gap 

W G (eV) 

0.1 

0.1 

0.1 

cr n (cm 2 ) 

10' 12 

io - 17 

2xl0' 16 

v p (cm 2 ) 

10 15 

10 13 

2xl0' 15 


Tableau IV.2 : Paramètres de la densité de défauts des matériaux en ZnO, CdS et CIGS 
[2-4], 


Après la détermination des modèles physiques vient l’étape de la définition de la méthode 
numérique de résolution et le calcul des paramètres de la cellule solaire. Nous avons choisi la 
méthode de Newton pour la résolution des équations de transport des semi-conducteurs liant le 
potentiel électrostatique et les densités de porteurs de charges ensemble. A chaque itération, le 
système des équations est résolu ensemble. La simulation électrique de la cellule solaire 
démarre par une polarisation nulle sur toutes les électrodes (état de l'équilibre 
thermodynamique). Une fois la solution de l'état de l'équilibre thermodynamique est trouvée, 
elle est utilisée pour l'initialisation de la simulation à l'état éclairée et polarisé. Les solutions 
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obtenues sont tracées soit par une interface graphique (TONYPLOT) d’Atlas ou par extraction 
des résultats. 

IV.3 Résultats de simulation de la cellule solaire en CdS/CIGS 
4- Cellule solaire en CdS/CIGS réalisée expérimentalement 

La cellule solaire en CIGS réalisée expérimentalement comporte sept couches minces. A 
l’arrière de la jonction, le substrat en verre (3 mm) et le molybdène (0.5-0. 9 pm) forme le 
contact ohmique arrière recueillant les trous. La troisième est la couche de CIGS (appelée 
absorbeur 2.5-3. 0 pm). La quatrième, appelée couche tampon, est une fine couche (30-50 nm) 
déposée à la surface du CIGS, et constituée de sulfure de cadmium CdS (de type n. La 
cinquième couche et la sixième couche sont à base de ZnO intrinsèque (50-100 nm) et ZnO 
dopé Al (150 - 200 nm) d’elle se situe à l’interface avec les couches suivantes à base de ZnO, 
puis de). La dernière couche en MgF 2 sert comme une couche anti réfléchissante de la lumière 
[1]. La figure IV.2(a) montre une coupe transversale et une coupe vue en microscopie 
électronique à balayage. 



Fig.IV.2 : Structure en coupe d’une cellule CIGS, (a) Image au microscope à Balayage, 
(b) schéma [12]. 
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La figure IV. 3 représente les caractéristiques I(V) et P(V) expérimentales d'une cellule 
solaire en CIGS avec une surface de 0.5 cm 2 . 
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Fig.IV.3 : Caractéristiques I(V) et P(V) expérimentales d'une cellule solaire en CIGS [1]. 


4 - Simulation de la cellule solaire en CdS/CIGS 

Après simulation de la cellule solaire en CdS/CIGS avec les paramètres précédents et une 
surface de 0.5 cm 2 en considérant le spectre solaire AM 1.5 de densité de puissance 100 
mW/cm 2 , nous avons obtenu les résultats représentés sur la figure IV.4. 
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Fig.IV.4 : Caractéristiques I(V) et P(V) simulées d'une cellule solaire en CIGS. 


Le tableau IV. 3 résume les paramètres photo voltaïques de la cellule solaire en CdS/CIGS 
simulés et ceux obtenus expérimentalement. 



/cc(mA. cm 2 ) 

V co ( mV) 

FF (%) 

7(%) 

Cellule réalisée 

expérimentalement 

37.8 

741 

80.6 

22.6 

Cellule simulé 

34.3 

803 

82.08 

22.6 


Tableau IV.3 : Comparaison entre les paramètres photo voltaïques simulés et 
obtenus expérimentalement des cellules solaires en CdS/ CIGS. 

On remarque que les résultats que nous avons obtenus par simulation sont en accord avec 
ceux trouvés expérimentalement. 

La cellule solaire en CdS/CIGS cause des problèmes dans l’environnement et la santé due à 
la toxicité du matériau en cadmium (Cd) contenu dans le sulfure de cadmium (CdS). On trouve 
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que les matériaux de sulfure de zinc (ZnS) sont les mieux convenables pour remplacer le sulfure 
de cadmium (CdS) [6,13-16], Ils sont dépourvu d’éléments toxiques et que la bande interdite 
du ZnS (3,68 eV) est supérieur à celle de CdS (2,42 eV), ce qui permet la transmission de 
photons à énergie supérieure à la couche absorbante en CIGS. C’est un matériau qui existe en 
abondance dans la nature et à faible coût, pour cela nous avons remplacé le matériau CdS par 
le matériau ZnS dans la suite de notre travail. 

IV.4 Résultats de simulation de Cellule solaire en ZnS/CIGS 

Nous avons pris une cellule solaire en ZnS/CIGS ayant la même structure que celle de la 
cellule solaire en CdS/CIGS décrite dans le paragraphe précédent (sous-titre IV.2.1), la couche 
en CdS a été substituée par la couche en ZnS. 

La figure IV.5 représente la structure simulée de la cellule solaire en ZnS/CIGS. 
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Fig.IV.5 : Structure bidimensionnelle d’une cellule solaire en ZnS/CIGS. 
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Les paramètres des matériaux ZnO, CdS, ZnS et CIGS utilisés dans la simulation de la 
cellule solaire sont reportés dans la littérature [1-6,17] et sont groupés dans le tableau IV.4. 


Régions 

Paramètres 

ZnO 

CdS 

ZnS 

CIGS 

Epaisseur (pm) 

0.05 

0.03 

0.03 

3 

E g (eV) 

3.3 

2.42 

3.68 

1.25 

N c (cm 3 ) 

2.2xl0 18 

2.2xl0 18 

1.5xl0 18 

2.2xl0 18 

N v (cm 3 ) 

1.8xl0 19 

1.8xl0 19 

1.8xl0 19 

1.8xl0 19 

Mobilité des électrons 
(cm 2 /V.s) 

100 

100 

250 

100 

Mobilité des trous 
(cm 2 /V.s) 

25 

25 

40 

25 

Permittivité relative 

9 

10 

8.32 

13.6 

Affinité 

électronique (eV) 

4.7 

4.5 

4.5 

4.8 


Tableau IV.4 : Paramètres de simulation la cellule solaire en CIGS. 

La cellule solaire en ZnS/CIGS sous éclairement par le spectre solaire AM 1.5 et de densité 
de puissance 100 mW/cm 2 et une surface de 0.5 cm 2 que nous avons simulé donne comme 
résultat : La densité du courant de court-circuit ) cc = 35.66 mA/cm 2 , la tension en circuit 
ouvert V co = 804 mV, le facteur de forme FF = 82.14 % et le rendement de conversion 
î] = 23.54 %. 
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La figure IV.6 représente les caractéristiques I(V) et P(V) d'une cellule solaire en ZnS/CIGS. 
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Fig.IV.6 : Caractéristiques I(V) et P(V) de simulation d’une cellule solaire en ZnS/CIGS. 

IV.5 Comparaison entre deux cellules en CdS/CIGS et ZnS/CIGS 

La figure IV.7 représente les caractéristiques J-V de deux cellules solaires en CIGS avec 
des couches tampon en CdS et en ZnS. 
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Fig.IV.7 : Courbe J-V de cellules en CdS/CIGS et en ZnS/CIGS. 

On remarque que la densité de courant de court-circuit (J cc ) de la cellule en ZnS est 
supérieure à celle de la cellule en CdS parce que l’absorption de la lumière pour les courtes 
longueurs d’ondes par le matériau CdS est supérieur à celle du matériau ZnS ( Ejj nS > E{j dS ). 
Par contre la tension en circuit ouvert ( V co ) de la cellule en CdS est presque la même que celle 
de la cellule en ZnS (Tableau IV.5). 



Jcc (mA. cm 2 ) 

Ko (y) 

FF(%) 

V (%) 

Cellule en CdS 

34.30 

0.803 

82.08 

22.60 

Cellule en ZnS 

35.66 

0.804 

82.14 

23.54 


Tableau IV.5 : Comparaison entre les cellules en CdS et en ZnS. 

La comparaison entre les deux cellules solaires en CdS et en ZnS indique que la cellule 
solaire en ZnS a de meilleure performance que celle en CdS. 
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Dans la suite nous allons étudier l'effet du dopage et de l’épaisseur des couches en ZnS et 
CIGS sur les paramètres de la cellule solaire en CIGS dans le but d’optimiser notre cellule 
solaire ZnS/CIGS. 

IV.5 Effets des couches en ZnS et CIGS sur les performances de la cellule solaire 
IV.5.1 Effet du dopage des couches en ZnS et CIGS 

Pour une couche en ZnS (30 nm) et CIGS (3 pm), nous avons calculé les paramètres de la 
cellule solaire pour différentes valeurs de la concentration du dopage N d du ZnS entre 
10 15 cm' 3 et lxlO 18 cm" 3 et avec différentes valeurs de la concentration du dopage N a du CIGS 
de 10 15 cm' 3 à lxlO 18 cm" 3 . Les résultats sont illustrés dans la figure IV. 8. 

Nous remarquons sur la figure IV. 8, que lorsque le dopage de la couche absorbeur en CIGS 
(. N a ) augmente de 10 15 à 10 18 cm" 3 , il y’a un décroissement de la densité du courant de court- 
circuit (J cc ) de 37.45 à 34.63 mA/cm 2 par contre il y a un faible effet des variations du dopage 
(j N d ) de la couche tampon en ZnS. Par contre la tension de circuit ouvert (V co ) augmente de 
682 à 882 mV quand le dopage ( N a ) de la couche absorbeur en CIGS augmente de 10 15 à 
10 18 cm" 3 alors que l’effet du dopage (N d ) de la couche tampon en ZnS est faible. 
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Fig.IV.8 : L’effet du dopage des couches en ZnS et CIGS sur les paramètres ( J cc , V co , FF, r\ ) de 
la cellule solaire. 
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Nous remarquons aussi une augmentation rapide du rendement Oj) pour un dopage entre 
2x1 0 17 et 10 18 cm' 3 des deux couches absorbeur en CIGS et tampon en ZnS causé par 
l’accroissement du champ électrique avec le dopage dans l’hétéroj onction ZnS/CIGS. Ceci 
permet une diminution de la recombinaison des photoporteurs dans la zone de charge d’espace, 
d’où une amélioration de la collecte des porteurs photogénérés et donc une augmentation du 
rendement de conversion. L’accroissement du rendement est principalement dû à 
l’augmentation de la tension en circuit ouvert et du facteur de forme (FF). 

Le rendement le plus élevé est de 26.06% obtenue pour une concentration de 6xl0 17 cm" 3 
de la couche en ZnS et une concentration du dopage de lxlO 18 cm' 3 de la couche en CIGS. 
Nous allons utiliser ces deux valeurs pour étudier l’effet de l’épaisseur des couches en ZnS et 
en CIGS sur les paramètres de la cellule solaire. 

IV.5.2 Effet de l’épaisseur des couches en ZnS et CIGS 

L'épaisseur de la couche en ZnS est variée de 10 nm à 50 nm tandis que l’épaisseur de la 
couche en CIGS est variée de 1 pm à 4 pm. L’effet de l’épaisseur des couches en ZnS et CIGS 
sur la cellule solaire est illustré dans la figure IV.9. 

Nous observons une augmentation de la densité da courant de court-circuit de 30.69 à 
35.70 mA/cm 2 avec l’accroissement de l’épaisseur de la couche absorbeur en CIGS de 1 à 
4 pm et pour les épaisseurs de la couche tampon en ZnS entre 10 et 35 nm. 
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Fig.IV.9 : L’ effet de l’épaisseur des couches en ZnS et CIGS sur les paramètres 
(J cc> Ko’ V ) de la cellule solaire. 


Chapitre IV 


Résultats et Discussions 


Nous remarquons aussi, que la tension en circuit ouvert V co augmente vers une valeur 
maximale V co = 0.882 V avec l’épaisseur de la couche tampon en ZnS de 0.03 pm et avec 
l’épaisseur de la couche absorbeur en CIGS de 4 pm. 

L’accroissement de V co et l cc engendrent une augmentation du rendement de conversion de 
la cellule solaire. Une couche absorbante très mince signifie physiquement que le contact en 
arrière et la région de déplétion sont très proches, ce qui favorise la capture d'électrons par ce 
contact. Cette forme de processus de recombinaison est préjudiciable aux performances de la 
cellule car il affecte V co , J cc et le rendement de conversion. 

La variation de l’épaisseur de la couche tampon entre 35 et 50 nm montre qu’il y a une 
diminution de tous les paramètres de la cellule (Jco^co’PP’V) avec l’augmentation de 
l’épaisseur de la couche ZnS. Le rendement maximal obtenu est 26.82 % avec une épaisseur 
de 25 nm de la couche tampon en ZnS et 4 pm de la couche absorbeur en CIGS. 

IV.6 Conclusion 

Dans ce travail, nous avons présenté et discuté les résultats d’étude par simulation 
numérique des caractéristiques électriques d’une cellule solaire en couches minces et à 
hétéroj onction à base de CIGS, ces résultats sont générés par le logiciel de simulation Silvaco 
Atlas-2D. 

On a comparé les résultats de la simulation électrique de la cellule en CdS avec ceux trouvés 
expérimentalement, ces résultats sont en accord avec ceux trouvés expérimentalement. En 
deuxième lieu, nous avons comparé les résultats de la simulation électrique de la cellule en 
CdS avec la cellule en ZnS, cette comparaison prouve que la cellule en ZnS est meilleure par 
rapport à la cellule en CdS. Ensuite nous avons étudié l’impact de deux couches en ZnS (couche 
tampon) et en CIGS (couche absorbeur) dans le but de concevoir une structure optimale à 
hétéroj onction ZnO/ZnS/CIGS qui donne le meilleur rendement électrique. 

Nous concluons de cette étude que le meilleur dopage de la couche en ZnS est de 
concentration 6x1 0 17 cm" 3 avec une épaisseur mince de 25 nm et le meilleur dopage pour la 
couche en CIGS de concentration lxlO 18 cm" 3 avec une épaisseur de 4 pm pour obtenir le 
rendement électrique optimal de 26.82%. 
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Conclusion générale 


Dans ce travail, nous avons étudié les caractéristiques électriques de la cellule solaire a 
hétéroj onction, ZnO/CdS/CIGS et ZnO/ZnS/CIGS, de structure (n-p). Ou la couche ZnO de 
type-n est utilisée comme une couche fenêtre, la couche CdS et ZnS de type-n comme une 
couche tampon et la couche CIGS de type-p est la couche de l’absorbeur, nous avons utilisé la 
simulation numérique par le logiciel Silvaco Atlas-2D, pour l’étude des caractéristiques de ces 
dispositifs. Nous avons aussi optimisé les paramètres de couches de la structure tels que 
l’épaisseur et la concentration du dopage pour obtenir un rendement de conversion électrique 
maximal. 

En première lieu, nous avons trouvé les résultats de simulation de la cellule solaire en 
ZnO/CdS/CIGS : La densité du courant de court-circuit ] cc = 34.30 mA/cm 2 , la tension en 
circuit ouvert V co = 0.803 V, le facteur de forme FF = 82.08 % et le rendement de 
conversion p = 22.60 %. Ces valeurs sont en accord avec celles trouvées expérimentalement. 

En deuxième lieu, nous avons comparé les résultats de simulation de la cellule solaire en 
ZnO/CdS/CIGS avec les résultats de la cellule solaire en ZnO/ZnS/CIGS, les résultats de cette 
dernière sont les suivants : La densité du courant de court-circuit ] cc = 35.66 mA/cm 2 , la 
tension en circuit ouvert V co = 0.804 V, le facteur de forme FF = 82.14 % et le rendement 
de conversion p = 23.54 %. Ces résultats prouvent que la cellule en ZnS est la meilleure. 

Par la suite nous avons étudié l’effet des deux couches ZnS et CIGS sur les performances 
de la cellule solaire (J cc , V co , FF, rj ). Nous sommes alors intéressés au dopage et à l’épaisseur 
des couches en ZnS et CIGS pour voir leurs influences sur le rendement électrique de la 
structure. 

En faisant varier le dopage la couche tampon ZnS de 10 15 à 10 18 et le dopage la couche 
absorbeur CIGS de 10 15 à 10 18 , un rendement optimal de 26.06% a été obtenu pour un dopage 
de 6xl0 17 cm' 3 pour ZnS et lxlO 18 cm' 3 pour CIGS. 

Et en faisant varier l’épaisseur de la couche tampon ZnS de 10 à 50 nm et l’épaisseur de la 
couche d’absorbeur CIGS de 1 à 4 pm, un rendement optimal de 26.82% a été obtenu pour une 
épaisseur de 25 nm pour ZnS et 4 pm pour CIGS. 
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En prenant une structure qui combine les paramètres optimaux (épaisseur et dopage) de 
chaque couche, nous avons pu atteindre un rendement de conversion électrique de 26.82% avec 
un facteur de forme de 85.37%. 

Les paramètres optimums de la structure optimale obtenus par notre simulation pour la 
couche tampon ZnS et la couche absorbeur CIGS sont présentés dans le tableau V suivant. 


ZnS 

CIGS 

n 

(%) 

FF 

(%) 

Ko 

(V) 

Jcc 

(mA/cm 2 ) 

N d 

( cm 3 ) 

Epaisseur 

(nm) 

N a 

i (cm 3 ) 

Epaisseur 

(pm) 

26.82 

85.37 

0.881 

35.64 

6xl0 17 

25 

IxlO 18 

4 


Tableau V : Paramètres optimum ayant donné r[=26.82 %. 

Cette étude nous a permis d’optimiser le rendement électrique de la structure 
ZnO/ZnS/CIGS et de montrer l’impact de la couche ZnS/CIGS sur les performances du 
dispositif en fonction du rôle de chacune d’elle. 
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